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1.  SUWARY 


Th«  objective  of  Che  TIES  (Tactical  Infocnatlon  Exchange  System)  program 
la  Che  Integration  of  Che  presently  used  separate  CNI  (Communication,  Naviga- 
tion and  Identification)  subsystems  into  a large  system  chat  is  progranaable, 
reconflgurable,  uses  comaon  modules,  and  can  process  a variety  of  signals 
simultaneously  that  are  currently  performed  one  at  a time  by  dedicated  hard- 
ware. The  system  Is  partitioned  into  three  parts  — Che  frequency  conversion 
unit,  the  frequency  distribution  unit,  and  the  signal  conversion  unit  as  shown 
In  figure  1. 

Because  of  the  difference  in  signal  structure,  the  signal  conversion  unit 
Is  further  broken  down  Into  the  narrow  band  unit  which  handles  narrow  band 
signals  such  as  AM,  IM,  SSB,  etc.,  and  the  wideband  unit  which  handles  signals 
such  as  TACAN,  IFF,  JTIDS,  and  GPS  In  Che  foreseeable  future. 

As  an  In  house  effort  Co  investigate  Che  narrowband  signal  processing  re- 
quirements for  TIES,  a narrowband  signal  conversion  unit  (NBSCU)  was  designed 
Chat  will  handle  AM,  FM,  and  SSB  signals  In  the  half  duplex  mode  using  presently 
available  digital  hardware.  The  unit  uses  FIR  (Finite  Impulse  Response)  filters 
because  FM  signals  are  extremely  phase  sensitive  and  FIR  filters  have  exactly 
linear  phase  response. 

The  Arctangent  algorithm  Is  used  for  signal  demodulation  because  of  Its 
efficiency.  The  amplitude  of  Che  signal  Is  given  by  Che  equation  A ■ Rcos6 
Qsln6.  The  angle  9 Is  given  by  6 ■ Cao'’^(Q/R).  In  both  cases,  R and  Q are 
Che  real  and  quadrature  components  of  Che  Incoming  signal  that  was  frequency 
translated  to  baseband. 

Detailed  descriptions  of  various  functions  to  be  performed  by  the  NBSCU; 
processing  algorlttiM  for  various  signals,  as  well  as  circuits  designed  to 
accomplish  these  functions  are  presented  In  Che  following  pages. 

2.  BASIC  CONCEPTS 
2.1  Sampling 

In  order  to  have  a better  understanding  of  the  work  Involved,  the  basic 
concept  of  sampling  has  to  be  reviewed.  The  sampling  theorem  states  that 
If  a band  limited  signal  Is  sampled  at  a sufficiently  high  rate,  Che  original 
signal  can  be  reconstructed  by  a set  of  Interpolating  functions  on  the  sampled 
values  of  the  signal.  The  mlninnim  sampling  rate  Is  the  well  known  Nyqulst 
race  of  twice  Che  highest  frequency  of  Che  Incoming  signal.  The  sampled  spec- 
trum of  a signal  Is  shown  In  figure  2. 

Undersampling  will  cause  the  repeated  spectrum  Co  fold  over  one  another 
causing  aliasing.  The  reconstruction  of  Che  original  signal,  depending  on  how 
severe  aliasing  Is,  becomes  extremely  difficult  and  can  only  be  recovered  with 
distortion.  The  amount  of  distortion  will  depend  on  Che  amotinC  of  under- 
sampling and  on  Che  shape  of  the  spectrum  of  the  signal  being  sampled.  This 
Is  shown  In  figure  3. 

If  Che  Incoming  signal  Is  at  a relatively  high  frequency,  sampling  at  twice 
Che  highest  frequency  will  become  impractical.  It  Is  known  Chat  for  a bandpass 
signal,  the  sampling  rate  can  be  as  small  as  twice  the  bandwidth  of  the  signal. 
Another  constraint  for  the  bandpass  sampling  frequency  Is  that  Che  sampling 
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Figure  2(a)  - Bandpass  Signal 

(b)  - Sampled  Signal  Spectrum 


Figture  3(a)  - Bandpass  Signal 

(b)  - Effect  of  Undersampling 
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race  muse  satisfy  Che  relationship 

F 

's  INT(Fh/B) 

where  is  che  highest  signal  frequency, 

B Is  che  bandwidth  of  che  signal, 

INT  Is  che  Integer  value  of  che  quanclCy. 

Therefore,  In  bandpass  sampling,  che  sampling  rate  varies  between  twice  Co  four 
Clmes  che  signal  bandwldch.  Another  effect  for  bandpass  sampling  Is  that  of 
frequency  translation.  This  Is  shown  In  figure  4. 

It  can  be  seen  then,  the  original  signal  can  be  recovered  independently  of 
che  sampling  technique  by  passing  the  sampled  data  through  a low  pass  filter 
with  che  proper  bandwldch  as  shown  In  figure  5. 

k 

It  Is  important  to  note  that  Che  duration  and  shape  of  che  pulses  in  * 

the  sampling  waveform  do  not  affect  the  ability  to  recover  the  signal  from  Its  , 

samples.  The  cask  for  processing  Che  conventional  signals  will  be  to  design 

a set  of  digital  hardware  that  will  perform  the  Interpolating  function  of  f 

input  samples  to  produce  che  original  baseband  signal.  The  state-of-the-art  t 

circuit  In  sampling  would  be  a sample-and-hold  circuit  which  consists  of  a high  | 

speed  electronic  switch  followed  by  a capacitor.  Upon  external  command,  Che  ! 

switch  Is  closed  and  the  capacitor  Is  charged  up  to  the  incoming  signal  level.  f 

The  switch  is  then  opened  and  Che  charge  is  then  held  on  by  che  capacitor  until  j 

che  next  consand  for  closing  the  switch  comes  along.  A simple  diagram  for  a 

sample-and-hold  circuit  is  shown  in  figure  6.  Sample-and-holds  normally  have 

unity  gain  and  are  non- Inver ting.  There  are  several  noise  sources  that  are  > 

Inherently  associated  with  sample-and-hold  circuits  that  will  be  discussed  in 

Appendix  A. 

2.2  A/D  Conversion 


The  sampled  output  of  che  signal  has  to  be  converted  into  a binary  number 
before  any  processing  can  be  performed.  This  process  is  known  as  A/D  conversion 
and  Is  accomplished  by  A/D  converters.  Most  commonly  used  A/D  converters  are  of 
2's  complement  notation  due  to  the  Inherent  advantage  of  the  2's  complement 
numbering  system  in  hardware  processing.  Since  che  input  voltage  to  the  con- 
verter needs  an  infinitely  large  number  of  binary  digits  for  representation, 
another  source  of  error  will  occur  due  to  the  finite  digit  representation  of 
che  analog  input  signal.  This  will  also  be  discussed  In  Appendix  A.  The 
output  digits  of  the  A/D  converter  will  be  the  inputs  to  the  narrowband  processor. 
To  reduce  circuit  complexity  to  a minimum  only  AM/FM/SSB  will  be  processed  one 
function  at  a time  In  the  half  duplex  mode. 

3.  SIGNAL  PROCESSING  ALGORITHMS 


3.1  AM  Transmit 


The  AM  transmit  algorithm  is  shown  in  figure  7. 
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Figure  4(a)  - Bandpass 


Figure  7 - AM  Transmit  Algorlttmi 
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Audio  Input  ranging  from  300  to  3300  Hz  Is  sampled  at  a rate  of  14.4KUz. 

The  sampled  Input  will  then  go  through  a digital  filter  that  passes  the  audio 
signal,  but  attenuates  signals  beyond  4KHz  to  50  db.  The  output  will  then  be 
Interpolated  by  two  Identical  1:4  digital  filters  to  change  the  sampling  rate  to 
230.4KHZ.  The  output  of  the  second  filter  will  be  at  audio  frequency  with  change 
In  sampling  rate  only.  A constant  number,  representing  the  carrier,  will  be 
added  to  this  audio  signal.  By  multiplying  the  resulting  samples  with  the 

sequence  1,  0,  -1,  0,  an  effective  mixing  process  with  samples  of  a cosine 

function  is  obtained.  The  mixing  process  causes  the  signal  to  shift  In  frequency 
by  ±57.6KHz,  corresponding  to  the  sum  and  difference  frequencies  in  analog 
mixing.  This  signal  la  then  transmitted  through  the  D/A  converter.  If  the 
carrier  Is  set  to  0,  then  double  sideband  suppressed  carrier  will  be  obtained. 

The  spectral  diagram  for  the  AM  transmit  chain  at  points  1,  2,  3,  4,  5,  6,  7 
Is  shown  in  figure  8.  From  the  diagram,  the  specifications  for  the  digital 
filters  can  be  obtained  and  are  shown  In  table  1. 


TABLE  1 - DIGITAL  FILTER  SPECIFICATIONS  FOR  AM  TRANSMIT  MODE 


Filter 

Sampling 

Frequency 

(KHz) 

Passband 

(KHz) 

Normalized 

Passband 

1 

14.4 

0-3.3 

0-0.2292 

2 

57.6 

0-3.3 

0-0.0573 

3 

230.4 

0-3.3 

0-0.01432 

3.2  AM 

Receive 

Stopband 

Normalized 

Attenuation 

(KHz) 

Stopband 

(db) 

4-7.2 

0.2780-0.5 

50 

10.4-28.8 

0.1806-0.5 

50 

53.6-115.2 

0.2326-0.5 

50 

The  AM  receive  algorithm  Is  shown  In  figure  9.  The  Input  signal  Is  con- 
verted Into  an  IF  frequency  of  57.6KHz  and  sampled  at  a rate  of  230.4KHz.  By 
multiplying  the  Incoming  signal  by  the  Incoming  IF  frequency  Is  then 

translated  to  baseband.  At  the  same  time,  two  channels  (one  real,  one  quad) 
are  formed.  This  can  be  seen  more  clearly  as  follows: 

Let  f(t)  be  the  Incoming  signal.  Multiplication  by  e'J^irfot  yields  the 
following: 


f(t)  cos  2TrfQt  - jf(t)  sin  2irfot 

Since  fg  Is  S7.6KHZ  while  the  Incoming  signal  Is  sampled  at  230.4KHz,  four 
samples  are  required  for  the  sine  and  cosine  functions  In  the  above  expression. 
For  ease  of  computation,  these  samples  can  be  chosen  to  have  the  following 
values : 

For  the  cosine  function  : 1,0, -1,0 

For  the  -sine  function  : 0,-1, 0,1 
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Therefore,  by  changing  Che  signs  of  Che  InpuC  samples  and  rouclng  Che 
samples  alcemacely  Co  Che  4:1  resampling  fllcer  wlch  approprlace  zero  filling 
as  shown  In  figure  9,  che  frequency  cranslaclon  can  be  accomplished  racher 
easily.  The  oucpuc  of  Che  frequency  cranslacer  will  feed  che  Cwo  4:1  de- 
sampling fllcers  on  che  real  and  quad  channels.  However,  close  examlnaclon  of 
che  block  diagram  Indlcaces  chac  che  funcclon  of  che  frequency  Cranslaclon  can 
be  readily  absorbed  InCo  chac  of  che  cwo  4:1  fllcers.  This  can  be  seen  from 
che  following: 


If  IC  Is  assumed  chaC  che  4:1  fllcer  Is  of  M scages  (where  N can  be  any 

Inceger)  wlch  coefflclenCS  for  che  Impulse  response  h^,  hj^ ^M-l* 

ouCpuC  of  che  fllcer  Co  a sec  of  InpuC  samples  Is  given  by  che  convoluClon  sum. 


N-1 

I Y(n)  - E h(k)  X (n-k) 

I k*o 

f 

%ihere  Y^Q^lsChe  ouCpuC  sample  and  X(n-k)  are  che  InpuC  samples.  This  process 
Is  shown  In  figure  10. 


However,  che  InpuC  daca  Co  Che  4:1  resampling  fllcer  In  chls  case  Is  mulclplled 

by  1,0, -1,0 for  che  real  channel  and  by  0,-1, 0,1 In  che  quadraCure 

channel.  Ic  can  be  seen  chaC  In  boch  channels,  alcemace  samples  have  value  0 
while  che  signs  change  alcemacely  for  samples  chac  are  noc  zero.  Using  chese 
InpuC  samples  for  Che  4:1  fllcer,  figure  11  Is  obcalned. 

Noclce  chac  InsCead  of  mulclplylng  Che  Input  samples  by  1,0,-1,0. . . .and  0, 
-1,0,1.... for  che  Cwo  channels,  che  coefflclenCs  for  Che  fllcer  can  be  varied 
Co  obcaln  ouCpuC  samples  as  shown  In  figure  12. 

Noclce  also  chac  while  Che  Cwo  ouCpuCs  of  che  cwo  figures  look  Idenclcal, 

Chey  are  In  face  noc  che  same  upon  close  examlnaclon.  A comparison  made  becween 
che  cwo  mechods  used  for  che  real  channel  shows  che  difference  for  Che  case  N>4. 
This  Is  shown  In  figure  13. 

Ic  can  be  seen  clearly  chaC  only  every  fourch  ouCpuC  of  che  Cwo  meChods 
are  Idenclcal.  While  che  ocher  cwo  secs  of  ouCpuCs  are  dlfferenc,  no  dlfflculcy 
Is  encouncered  because  In  a 4:1  fllcerlng  process  only  one  ouc  of  every  fourch 
consecuClve  ouCpuC  Is  used.  Thus,  for  chls  appllcaClon,  che  funcclons  for  che 
Cwo  mechods  are  Idenclcal.  Consequencly,  che  frequency  cranslaclon  hardware 
can  be  ellmlnaCed.  The  fllcerlng  process  can  be  simplified  by  combining  che 
Cwo  fllcers  In  che  Cwo  channels  InCo  a single  fllcer.  This  can  be  accomplished 
by  changing  che  sign  of  alcemace  pairs  of  coefflclenCs  and  rouclng  alcemace 
sum  of  produces  Co  cwo  secs  of  reglscers  for  s Corage.  This  Is  shown  In  figure  14. 
While  chls  looks  more  compllcaced  chan  che  single  channel  operaclon,  Ic  Is  In 
face  very  simple  Co  ImplemenC  as  will  be  seen  lacer  on. 

The  baseband  signal,  due  Co  che  process  of  frequency  cranslaclon,  will  be 
complex.  Ic  has  a real  componenc  as  well  as  a quadraCure  componenC  as  shown 
In  figure  IS. 

The  magnlCude  of  che  veccor  will  Chen  be  che  desired  AH  oucpuc.  Two  mechods 
can  be  used  for  chls  compuCaClon  and  are  given  by  che  following  Cwo  equaciona. 
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(a)  Real  Channel 

(b)  Quad  Channel 
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where  A - Amplitude  of  vector 

R ■ Real  component  of  vector 
Q > Quadrature  component  of  vector 
9 > Phase  angle  of  vector 


I 

f 


The  first  equation  is  obtained  from  Pythagoras  Theorem  while  the  second 
equation  is  obtained  from  the  trigonometric  relationship  between  the  sine  and 
cosine  functions  as  follows: 


cos  9 ■ R/A 
sin  9 • Q/A 
Therefore, 

co8*9  ■ R cos  9 /A 
sln^9  ■ Q sin  9 /A 

Adding  Che  two  preceding  equations  yields: 
1 ■ (R  cos  9 + Q sin  9 ) /A 


Multiplication  by  A on  both  sides  of  Che  expression  yields  the  following 
equation: 

A > R cos  9 + Q sin  9 

Since  taking  the  square  root  of  a number  is  a complicated  process,  the 
second  method  is  chosen.  Notice  that  A is  Che  amplitude  of  a vector  and  therefore, 
the  sign  is  always  positive.  The  spectra  at  points  1,2, 3, A, 5 for  the  AM  receive 
mode  is  shown  in  figure  16.  The  specifications  for  Che  4:1  resampling  filter 
are  listed  in  cable  2. 


TABLE  2 - DIGITAL  FILTER  SPECIFICATIONS  FOR  AM  RECEIVE  MODE 
Sampling 

Frequency  Passband  Normalized  Stopband  Normalized  Attenuation 
Filter  (KHz)  (KHz)  Passband  (KHz)  Stopband  (db) 


230.4 


0-3.3  0-0.01432  53.6-115.2  0.2326-0.5 
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3.3  FM  Transmit 


FM  Is  a form  of  exponential  modulation.  In  exponential  modulation,  the 
modulated  wave  In  phasor  form  Is  an  exponential  function  of  the  message;  that 
is, 

Xg(t)  ■ Re  [Acejec(t)] 

- Ac  cos  0c(t) 

where  6c(t)  is  a linear  function  of  x(t),  the  baseband  message,  and  A^  is 
constant.  Furthermore, 

ec(t)  - Zirfgt  + (j»(t), 

so  that  the  carrier  frequency  f^  Is  uniquely  specified.  The  second  term  in  the 
equation,  (^(t),  will  then  be  the  relative  phase  angle  In  the  sense  that  the 
phasor  ej^c  differs  In  angular  position  from  eJ^ct  by  ^(t).  Since  frequency  Is 
the  time  derivative  of  angular  position,  the  instantaneous  frequency  deviation 
f(t)  can  be  written  as 


f(t) 


1 d»(t) 
iTT  dt 


Therefore , 

ec(t)  - 2Trfgt  + 2ir  /fCDdX 


But  the  instantaneous  frequency  Is  also  proportional  to  the  message. 

Therefore, 

f(t)  - f^X(t) 

The  equation  for  6c(c)  can  be  rewritten  as  follows: 

ec(t)  - 2irfc(t)  + 2Trf^/xV)dX 

*00 

The  modulated  wave  will  then  be: 

Xc(t)  ■ Re  [AceJ®c(t)] 

- Re  Ac  [cos  0c(t)  + j sin  0c(t)] 

■ Ac  cos[<0ct  + 2ltf^  j x(X)dA] 

*00 

The  algorithm  for  FM  transmit  Is  shown  In  figure  17. 

Audio  Input  is  sampled  at  a rate  of  14.4KHz  and  the  sampling  rate  Is  up 
converted  to  230.4KHz  after  filtering  and  resampling.  The  output  of  the  filter 
Is  then  Integrated  to  obtain  the  phase.  The  Integration  process  Is  given  by 
the  following  equation: 

Y(n)  - X(n)  + Y(n-l) 
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Where  Y(n)  Is  the  present  output  sample,  X(n)  Is  the  present  Input  sample,  and 
Y(n-l)  Is  the  delayed  output  sample. 

Recall  that 

0c(t)  - 2iifc(t)  + 2:Tf^  /xCDdA 

im 

The  product  2iTf^  Is  a scaler  and  can  be  Incorporated  into  the  filter  coefficients. 

In  order  to  obtain  the  desired  FM  signal,  the  real  part  has  to  be  extracted 
from  the  function  as  follows: 

FM  - Re  [eJ®c(0] 

• Re  j x(X)dXj 

Let  f x(X)dA  . gj<|, 

■ cos  ♦ + j sin 

and  eJ2xfc(t)  . 2Trfc(t)  + J sin  2Trfc(t) 

Therefore,  ^ 

Re  [eJ2fffc(t)gJ2irf/i  J x(A)dA] 

• Re  [cos  2irfg(t)  + J sin  27rfg(t)][co8  ♦ + J sin  ♦] 

- cos  4)  cos  2Trfg(t)  - sin  (>  sin  2irfg(t) 

If  the  carrier  frequency  is  chosen  to  be  a quarter  of  the  input  sampling 
tate,  then,  four  samples  of  the  cosine  function  (which  can  be  chosen  to  be 
1,0,-1,0, . . . . ) and  four  samples  of  the  sine  function  (which  can  be  chosen  to  be 

0«1>0,-1, ) can  be  used  for  multiplication  to  produce  the  desired  output 

as  follows: 

(co8<^) [l,0,— 1,0, . . . . . ^ — (sin  4^)  [0,1,0,— 1, . . . . 3 

■ cos  (p,  - sin  (p,  - cos  (p,  + sin  (p,  + 

Therefore,  the  output,  Y(n)  (corresponding  to  ^ in  the  above  expression), 
from  the  Integration  process  can  be  used  as  an  address  to  obtain  the  sine  and 
cosine  values  of  Y(n).  By  changing  the  sign  of  every  other  pair  of  samples, 
the  frequency  translation  is  completed.  The  FM  transmit  spectra  is  shown  in 
figure  18.  The  specifications  for  the  filters  are  listed  in  table  3. 
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Figure  18  - FM  Transmit  Spectra 
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TABLE  3 - DIGITAL  FILTER  SPECIFICATIONS  FOR  FM  TRANSMIT  MODE 

Sanpllng 


Filter 

Frequency 

(KHz) 

Passband 

(KHz) 

Normalized 

Passband 

S topband 
(KHz) 

Normalized 

Stopband 

Attenuation 

(db) 

1 

14.4 

0-3.3 

0-0.2292 

4-7.2 

0.2780-0.5 

SO 

2 

57.6 

0-3.3 

0-0.0573 

10.4-28.8 

0.1806-0.5 

50 

3 

230.4 

0-3.3 

0-0.0143 

53.6-115.2 

0.2326-0.5 

50 

3.4  FM  Receive 

The  algorithm  for  the  FM  receive  mode  is  shown  in  figure  19.  As  can  be 
seen,  it  Is  Identical  to  the  AM  receive  mode  with  the  addition  of  a differen- 
tiator at  the  output  of  the  arctangent  demodulator.  The  Input  to  the  differ- 
entiator Is  the  phase  output  of  the  arctangent  demodulator.  Output  of  the 
differentiator  will  be  that  of  frequency.  The  FM  receive  spectra  at  points 
1,2, 3, 4, 5, 6 are  shown  In  the  FM  receive  algorithm  In  figure  20. 

The  specifications  for  the  filters  will  be  the  same  as  those  for  the  AM 
receive  mode.  As  for  the  differentiator.  In  the  frequency  domain,  it  will  have 
a linear  slope  in  the  frequency  band  of  Interest. 

3.5  Single  Side  Band  Transmit 

The  single  sideband  signal  can  be  obtained  from  the  double  sideband  signal 
by  suppressing  the  carrier  and  filtering  out  one  of  the  sidebands  (upper  or 
lower  sideband).  The  suppression  of  the  carrier  can  be  easily  done  by  adjusting 
the  constant  added  to  the  output  waveform  before  frequency  translation.  How- 
ever, since  the  two  sidebands  are  only  several  hundred  hertz  apart,  the  filter 
has  to  be  extremely  selective  and  represents  a relatively  large  load  on  the 
multipliers.  Consequently,  this  approach  is  abandoned. 

The  algorithm  for  single  sideband  transmit  Is  shown  on  figure  21.  Audio 
signals  ranging  from  300  to  3300  Hz  are  sampled  at  7.2KHz  and  frequency 
translated  1800  Hz  to  DC  frequency  by  multiplying  the  input  by  eij2Tr(1800) . 

The  positive  sign  Indicates  a frequency  up  shift  to  the  right  of  the  DC 
frequency  while  a negative  sign  indicates  a down  shift  to  the  left  of  the 
frequency  origin.  Due  to  the  frequency  translation  process,  the  signal  Is 
split  into  the  real  and  quadrature  components.  A low  pass  selectivity  filter 
will  then  eliminate  one  of  the  sidebands  depending  on  vdilch  way  the  frequency 
translation  took  place.  This  can  be  seen  more  clearly  on  the  spectra  diagram 
for  single  sideband  transmit  on  figure  22  in  which  the  spectra  at  points  1,2, 

3, 4, 5, 6, 7, 8 on  the  SSB  transmit  algorithm  are  shown. 

The  signal  then  goes  through  three  interpolating  filters  to  convert  the 
sampling  rate  to  230.4KBz.  By  multiplying  the  signal  with  4 samples  of 
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Figure  22  - SSB  Transmit  Spectra 
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e~^  t a frequency  translation  of  57.6  KHz  Is  accomplished.  The  real 

part  of  this  signal  is  the  desired  single  sideband  waveform.  To  determine 
whether  upper  sideband  or  lower  sideband  is  used,  one  of  the  frequency  trans- 
lation processes  is  chosen  to  have  a fixed  shift  (up  or  down).  If  it  is 
assumed  that  the  57.6  KHz  translation  has  an  up  shift,  a down  shift  for  the 
1.8  KHz  translation  will  produce  upper  sideband  while  an  up  shift  of  1.8  KHz 
will  produce  lower  sideband.  Again,  the  samples  for  frequency  translation 

are  chosen  to  be  the  two  sets  of  samples  (1,0, -1,0 0,1,0,-!,....)  for 

the  sine  and  cosine  functions  so  that  multiplication  is  reduced  to  a minimum. 

For  the  1.8  KHz  translation,  we  have  the  following: 

X(n)[cos  2nl800t  S j sin  2xl800t] 

- X(n)[(l,0,-1,0 ) i j (0,1, 0,-1, )] 

which  is  complex. 

The  case  for  the  57.6  KHz  translation  is  different  because  the  input  to 
the  frequency  translator  is  a complex  number.  Let  the  input  complex  nus^ers 
be  represented  by  x ■f  Jy.  The  output  of  the  frequency  translator  is  then  given 
by  the  following: 

(*  Jy)[cos  (2x57600) t + j sin  (2x57600) t] 

Taking  the  real  parts  yield 

X cos  (2x57600) t - y sin  (2x57600) t 

Therefore,  only  every  other  sample  in  each  channel  Is  needed  to  produce  the 
output  if  four  samples  of  sine  and  cosine  are  used. 

Notice  that  the  final  output  is  not  the  same  as  that  of  AM  double  side- 
band with  the  carrier  and  one  sideband  filtered  out.  This  can  be  seen  clearly 
in  figure  23.  The  difference  being  that,  in  using  this  method  of  modulation, 
the  transmitted  signal  is  centered  on  57.6  KHz,  while  filtering  the  carrier 
and  one  sideband  of  the  double  sideband  AM  signal  produce  a signal  that  is 
offset  from  the  carrier  frequency  of  57.6  KHz.  If  a cone  is  used  as  an  example, 
AM  modulation  will  produce  the  spectrum  shown  in  figure  23(a).  SSB  modulation 
which  will  be  used  for  the  MBSCU  will  produce  the  spectrum  in  figure  23(b). 

By  filtering  the  AM  signal,  Che  resulting  spectrum  is  shown  in  figure  23(c). 

The  difference  is  obvious.  The  filter  specifications  are  listed  in  Cable  A. 


TABLE  4 - DIGITAL  FILTER  SPECIFICATIONS  FOR  SSB  TRANSMIT  MODE 


Sampling 

Frequency  Passband  Normalized  Stopband  Normalized  Attenuation 
Filter  (IMz)  (KHz)  Passband  (KHz)  Stopband  (db) 


1 

7.2 

0-1.65 

0-0.2292 

2 

14.4 

0-1.65 

0-0.1146 

3 

57.6 

0-1.65 

0-0.0287 

4 

230.4 

0-1.65 

0-0.0072 

2-3.6  0.2780-0.5  50 

5. 2-7. 2 0.3610-0.5  50 

12.4-28.8  0.2153-0.5  50 

55.6-115.2  0.2410-0.5  50 


i 
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Figure  23  - (a)  DSB  Spectra 

(b)  SSB  Spectra 

(c)  SSB  Spectra  Obtained  by  Filtering  DSB  Signal 
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3.6  Single  Sideband  Receive 

The  single  sideband  receive  mode  Is  processed  similarly  Co  that  of  single 
sideband  cransmlc  mode  except  In  Che  reverse  order.  The  algorithm  Is  shown  In 
figure  24  and  Che  spectra  at  points  1,2, 3, 4, 3,6, 7 are  shown  In  figure  25. 

If  the  57.6  KHz  frequency  translation  Is  an  up  shift  operation,  the  result 
would  be  lower  sideband.  If  It  Is  a down  shift  operation,  Che  result  would  be 
upper  sideband  provided  Chat  the  1.8  KHz  frequency  translation  is  an  up  shift 
operation  In  both  cases.  The  filter  specifications  are  shown  In  cable  5. 

TABLE  5 - DIGITAL  FILTER  SPECIFICATIONS  FOR  SSB  RECEIVE  MODE 
Sampling 


Filter 

Frequency 

(KHz) 

Passband 

(KHz) 

Normalized 

Passband 

Stopband 

(KHz) 

Normalized 

Stopband 

Attenuation 

(db) 

1 

230.4 

0-1.65 

0-0.0072 

55.6-115.2 

0.2410-0.5 

50 

2 

57.6 

0-1.65 

0-0.0287 

12.4-28.8 

0.2153-015 

50 

3 

14.4 

0-1.65 

0-0.1146 

5. 2-3. 6 

0.3610-0.5 

50 

4.  CRITICAL  FUNCTION  DESIGNS 

As  can  be  seen  from  all  Che  algorithms,  the  most  Important  tasks  the 
NBSCU  has  to  perform  are  those  of  filtering  and  arctangent  demodulation.  These 
are  discussed  In  the  following  paragraphs. 

4.1  Digital  Filter  Design 

The  design  of  the  digital  filter  Is  Che  most  Important  Cask  for  the 
NBSCU  as  It  Is  used  In  all  the  mod/demod  algorithms.  In  general,  a digital 
filter  performs  the  equivalent  of  an  analog  filter  In  that  It  Is  frequency 
selective.  However,  the  same  digital  filter  can  be  made  to  be  a differentiator 
or  a Hilbert  transformer  by  varying  a set  of  tap  weights  for  the  filter. 

Unlike  analog  filters,  which  can  be  designed  to  pass  certain  frequency  bands 
from  DC  to  Infinity,  digital  filters  do  not  enjoy  this  luxury.  Ml  frequencies 
are  normalized  Co  Che  sampling  frequency  F^,  which  Is  taken  to  be  1 (or  ir)  In 
Che  digital  domain  and  frequencies  normalized  to  F^  are  known  as  digital  fre- 
quencies. Due  to  spectral  folding  In  a sampled  system,  the  Input  Co  a digital 
filter  must  not  have  any  significant  spectral  components  beyond  F3/2.  Other- 
wise, Che  original  signal  cannot  be  recovered  without  distortion  (the  sampling 
theorem).  Therefore,  the  analog  Input  to  an  A/D  converter  must  be  band  limited 
by  an  anti-aliasing  filter. 

There  are.  In  general,  two  classes  of  digital  filters: 

1.  Infinite  Impulse  response  (HR)  filters  in  which  the  impulse  response 
of  Che  filter  Is  Infinitely  long,  and 


Figure  24  - SSB  Receive  Algorithm 


Riacelve  Spectra 
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2.  Finite  Impulse  response  (FIR)  filters,  which  have  finite  time  Impulse 
responses  and  can  be  obtained  by  truncating  the  Infinite  Impulse  response 
samples. 

Analog  filters  such  as  Butterworth,  Chebyshev,  and  Elliptic  filters  which 
have  polynomials  both  In  the  denominator  and  numerator  In  the  transfer  function 
In  the  S plane  belong  to  the  first  class.  Their  digital  filter  counterparts 
which  have  the  same  kind  of  frequency  response,  can  be  obtained  by  using  a 
bilinear  transformation.  The  results  are  excellent  using  this  method.  How- 
ever, the  Implementation  Is  recursive,  l.e.,  the  output  samples  are  dependent  I 

on  the  present  Input  samples  as  well  as  the  previous  output  samples  and  care 
must  be  taken  that  limit  cycles  do  not  occur*.  Phasewlse,  the  HR  filter 
response  Is  extremely  nonlinear.  A second  order  section  of  an  HR  filter  Is 
shown  In  figure  26.  The  letters  a,  b,  c,  d.  In  the  figure  are  the  coefficients  i 

for  the  second  order  polynomial. 

b 

To  alleviate  the  problem  of  limit  cycles  and  phase  non-llnearltles,  the  . 

FIR  approach  can  be  taken.  To  generate  a FIR  filter,  the  Impulse  response 
of  a certain  filter  with  the  desired  characteristic  Is  obtained.  (These  may 
be  the  truncated  impulse  responses  from  HR  filters).  A convolution  between  r' 

the  Impulse  response  of  the  filter  and  the  Input  samples  to  the  filter  will  ^ 

produce  the  desired  output.  FIR  filters  have  exactly  linear  phase  and  since  |[; 

; a convolution  Is  a sum  of  product  computation,  the  absence  of  Input  will  pro-  f 

duce  no  output.  Hence,  there  Is  no  limit  cycle.  There  are  many  ways  In  which  t 

FIR  filters  can  be  designed  (besides  truncating  the  HR  filter  Impulse  response)  ^ 

such  as  the  use  of  various  windows,  the  frequency  sampling  technique ^ etc. 

A design  technique  which  Is  widely  accepted  In  recent  years  to  yleli  an 
optimum  (mlnlmax  error)  FIR  filter  Is  known  as  the  Remez  Exchange  algorithm. 

It  proved  to  be  the  best  of  the  various  design  techniques  listed  In  Appendix  B. 

In  general,  HR  filters  are  more  efficient  than  the  FIR  filters  when  the 
order  of  the  filter  Is  large.  However,  since  FIR  filters  have  exactly  linear 
phase,  l.e.,  there  Is  no  group  delay  distortion,  FIR  filters  are  chosen  for 
the  NBSCU. 

The  Remez  Exchange  algorithm  Is  Implemented  using  the  CDC  6600  to  obtain 
the  filter  coefficients  which  are  listed  In  tables  6 through  11. 

The  frequency  response  using  these  coefficients  are  then  obtained  using 
an  Interpolation  formula  Implemented  on  a Tektronix  4051  terminal.  The  output 
samples  are  plotted  and  shown  In  figure  27  through  figure  32  (with  2 db/dlvlslon 
vertical  scale).  The  response  with  the  coefficients  rounded  off  to  12  bits 
are  shown  In  figure  33  through  figure  38.  Notice  that  In  all  cases,  the 
attenuation  In  the  stopband  Is  at  least  50  db  down  with  respect  to  passband 
^ amplitude. 


*Llmlt  cycle  Is  an  effect  due  to  the  dependence  on  previous  output  samples. 
The  filter  produces  an  output  sample  with  no  Input. 
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TABL£  6 


- IMPULSE  RESPONSE  FOR  10  STAGE  1:2  DIGITAL  FILTER 


FINITC  IMPULSE  PESPONSS  IFIR) 

- LINEAR  PHASE.  31  GITAL -FILTER-  OESIGN 

REMcZ  EXCHANGE  'ALGORITHM 

--  aANLP ASS- filter  _ 


filter  length  = 10 

FILTER  LENGTH  CcTERMINED  BY  APPROXIMATION 


lower  banc  edge 

UPPER  BAND  EDGE 
-CESIREC  VALUE 
WEIGHTING 
DEVIATION 
QE-VIATICN-IN-Qi- 


XMEULSE-RESRONSE- 


H( 

1) 

-.13634021E-C1 

s 

HI 

10) 

H( 

2) 

5 

“ . 40  3 680  55  E-P 1 

HI 

9) 

HI 

3) 

S 

-.14243194E-02 

= 

H 1 

6) 

HI 

41 

r 

.I7ei8923£^00 

s 

HI 

71 

HI 

5» 

s 

.37921318E+00 

= 

HI 

6) 

■ -BANO 

-1 BANC  . 

-2- 

.OOOOOCOOO  .361300000 

.1146000C0  .50C000C00 

— 1^004-0044100 .4000404100 


1.0  000C300C  lOO.OOC  OOJCJ'i 
.0287<3210C  .300287931 

->5403256 49 -7 (UOIL 43LS5L  


band 


EXTREMAL  FREQUENCIES 

>0642540 >1.146000 

.4735000 


42 
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TABLE  7 - IMPULSE  RESPONSE  FOR  10  STAGE  1:4  DIGITAL  FILTER 


FINITE  IMPULSE  RESPONSE  (FIR) 

-LXNSAiL  -FHASE-DIilXI  Al,*  -FXL-TiR  --DESIGN 

REMEZ  EXCHANGE  ALGORITHf' 

-S,ANCPASS-  -FILTER 

filter  LENGTH  = IF 

FILTER  LENGTH  DETERMINED  BY  APPROX  I hATICN 


LOWER  BANC  EOGE 
UPPER  BAND  EDGE 


WEIGHTING 
DEVIATION 
-CEV. 


1) 

s 

meUCSE  RESPONSE 
.823346A1E-02 

= 

H ( 

1C) 

H( 

2) 

2 

.34954g63E-01 

2 

H( 

9) 

H( 

3) 

S 

.84412599E-C1 

3 

H( 

8) 

HC 

4) 

= 

.14330758E+00 

S 

H( 

7) 

K 

5) 

s 

.184197e4F^0F 

S 

H ( 

6) 

-flAt 

ID 

J aAND 

_2_ 

.000303000 

.007200300 


1 1 i ' 11 II 1 1 Irll 1 1 1 (■ 


i.COOOCQJOO 

.CS2834110 


.241000000 

.500000000 


IrVtiiWiiMiTil 


ICC.OOCOOCCOO 

.000928341 


EXTREMAL  FREQUENCIES 

.UlflZXGGO .^411 

.4597500 


— 84-7  500 


I 
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TABLE  8 - IMPULSE  RESPONSE  FOR  25  STAGE  1:4  DIGITAL  FILTER 


FINITE  IMPULSE  RESPONSE  CFIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 


BANW»ASS  filter 


FILTER  LENGTH  = 25 

IMPULSE  RESPONSE  ..... 

H< 

1)  = 

. 2110C766E-C3  = HI 

25) 

H( 

2)  = 

-.390C9108E-03  = HI 

24) 

HI 

3)  = 

-.28928785E-C2  = HI 

23) 

HI 

4)  = 

-. 8359993 8E-02  = HI 

22) 

HI 

5)  = 

-.16C27747E-C1  = HI 

21) 

HI 

6)  = 

-.22021152E-01  = HI 

20 

HI 

7)  = 

-.1954131 4E-01  = HI 

19) 

HI 

81  = 

-.1I701737E-02  = HI 

181 

HI 

91  = 

.35973699E-01  = HI 

17) 

HI 

10)  = 

.88589167E-01  = HI 

16) 

HI 

11)  = 

.1445e815E+00  = HI 

15) 

HI 

12)  = 

.18771701E+00  = HI 

14) 

HI 

13)  = 

.2039388CE4-00  = HI 

13) 

BAND 

1 BAND  2 

BAND 

LONER  BAND 

EDGE 

.000008000  .180600000 

UPPER  BAND 

EDGE 

.057300000  .500000000 

DESIRED  VALUE 

1 

.000000000  .000000000 

WE  IGHTING 

1 

.000000000  100.000000000 

DEVIATION 

.023109854  .000231099 

DEVIATION 

IN  OB 

.401530770  -72.724056063 

EXTREMAL  FREQUENCIES 

.OOOCOC3 

.0  384615 

.0573000  .1806000 

.1 878115 

.2070423  .2358885  .2695423  .3056000  .3440615 

._3825231  _ .4  20  9 8 46  . ..4  5 944  6 2 . 5 0 0 0 0CO  _ _ 


r 


1-1 
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TABLE  9 - IMPULSE  RESPONSE  FOR  55  STAGE  1:1  DIGITAL  FILTER 

FINITE  IMPULSE  RESPONSE  (FIR) 

LINEAR  PHASE  DIGITAL  FILTER  DESIGN 

REMEZ  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 

FILTER  LENGTH  = 55 


IMPULSE  RESPONSE 


H( 

1)  = 

.96952259E-03 

= 

HC 

55) 

H( 

2)  = 

. 10626658E-02 

S 

HC 

54) 

H( 

3)  = 

-.2798C042E-C2 

= 

H C 

53) 

H( 

4)  = 

80329778E-C2 

s 

H ( 

52) 

Hf 

5)  = 

-.6297514CE-C2 

= 

HC 

51) 

H( 

6)  = 

.22990183E-02 

= 

HC 

50) 

H( 

7)  = 

.47051969E-C2 

s 

HC 

49) 

H( 

8)  = 

-.36383487E-C2 

s 

H ( 

48) 

H( 

9)  = 

-.67g2C77QE-C2 

r 

H C 

47) 

H( 

10)  = 

.39219972E-02 

= 

HC 

46) 

H( 

11)  = 

. 9C7C4774E-02 

HC 

45) 

H( 

12)  = 

-.44368474E-fl2 

= 

H ( 

44) 

H( 

13)  = 

-.12204033E-C1 

3 

HC 

43) 

H( 

14)  = 

.48794071E-02 

=: 

HC 

42) 

H( 

15)  = 

.16271823E-01 

3 

HC 

41) 

H( 

16)  = 

52916459E-C2 

s 

HC 

40) 

HI 

173  = 

-.21694685E-01 

3 

HC 

39) 

H( 

18)  = 

•56489810E~02 

3 

HC 

38) 

H« 

19)  = 

.29229548E-01 

3 

HC 

37) 

HC 

20)  = 

-.59589844E-C2 

= 

HC 

36) 

H{ 

21)  = 

-.4(I564246E-01 

3 

HC 

35) 

H( 

22)  = 

.62099160E-02 

3 

HC 

34) 

HC 

23)  = 

.60088016E-01 

= 

HC 

33) 

H( 

2 4)  = 

-.638448C5E-C2 

= 

H( 

32) 

H( 

25)  s 

-.10392688E+00 

3 

HC 

31) 

H( 

26)  = 

.64950402E-02 

= 

HC 

30) 

HC 

27)  = 

.31757675E+CC 

3 

HC 

29) 

HC 

2 8)  = 

.49346309E+C0 

= 

HC 

28) 

BAND 

1 BAND 

2 

BA 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 

.000000000  .278000000 

.22923COOO  .503000000 

l.CCOCOOOOO  .000000000 

WEIGHTING  l.CCOncCOOO  100.000000030 
DEVIATION  .025721650  .000257216 
DEVIATION  IN  DP  .446929404  -71.794023574 

EXTREMAL  FREQUENCIES 

.CCOOOOJ 

.0212054 

.0424107  .0636161 

.0848214 

.104  9107 

.1250000 

.1462054  .1662946 

.185  26  79 

.2053571 

.2209821 

.2292000  .2780000 

.2813482 

.2891607 

.3C14375 

.3148304  .3304554 

.346C8C4 

.361  7C54 
.4476429 

.3784464 

.4655000 

.39630  36  .4130446 
.4822411  .5000000 
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TABLE  10  - IMPULSE  RESPONSE  FOR  16  STAGE  1:4  (SSB)  DIGITAL  FILTER 


PINITE  IMPULSE  RESPONSE  (PIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REME2  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 


FILTER  LENGTH  = 16 


IMPULSE  RESPONSE 


LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
OESIREO_yALUE 
HEIGHTEN G 
DEVIATION 
DEVIATION  IN  DB 


H( 

11  = 

-.92734092E-C3 

H { 

16) 

'H{' 

2)  = 

-.12691651E-02 

= 

H( 

15) 

HC 

31  = 

. 33726088E-C2 

= 

H ( 

14) 

HC 

41  = 

.20226966E-01 

s 

HI 

13) 

H{ 

5)  =“ 

.54508616E-C1 

s 

HI 

12) 

H( 

6)  = 

. 10325178E+00 

s 

H ( 

11) 

H( 

71  = 

.15292869E+00 

s 

HI 

1C) 

HI 

81  = 

.18463465E*^00 

s 

HI 

9) 

BAND 

1 BAND 

2 

.EQOCCCOOO 

.026700000 

1^.000000000 

l.COOOOCflOO' 

.033453619 

.581365829 


.2153COCOO 

.500000000 

.OCCOOOQOO 

100.000000030 

.000334536 

-69.511137847 


BAND 


EXTREMAL  FREQUENCIES 

_ ,.0Ca00C3_  .0  287000  __  .2153000  .2270188 

.3051437  .3598312  .4145187  ' .4692062 


.2621750 
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TABLE  11  - IMPULSE  RESPONSE  FOR  16  STAGE  1:4  DIGITAL  FILTER 


FINITE  IMPULSE  RESPONSE  (FIR) 

LINEAR  PHASE  DIGITAL  FILTFR  DESIGN 

REMEZ  EXCHANGE  ALGORITHM 

BANDPASS  filter 

FILTER  LENGTH  = 16 

• 

H( 

1)  = -.27836378E-03  = H( 

16) 

H( 

2)  = .19176790E-03  = H( 

15) 

H( 

3)  = .5138C971E-02  = H( 

14) 

H( 

4»  = .21038964E-01  = H( 

13) 

HC 

5»  = .53192579E-01  = H( 

12) 

HI 

6)  = . 99517446E-C1  = H( 

11) 

H( 

7)  = .14739916E«-00  = H( 

10) 

H( 

8)  = .17825397E+00  = H( 

9) 

BAND  1 BAND  2 

BAND 

LOWER  BAND  EDGE 

.CCOCOOOOO  .232600000 

UPPER  BAND  EDGE 

.014300000  .500000000 

DESIRED  VALUE 

1 

.000000000  .000000000 

WEIGHTING 

1 

.CCOCOOOOO  lOO.OCCOOOGOG 

DEVIATION 

.008907236  .000089072 

DEVIATION  IN  DB 

.154738637  -81.005140433 

EXTREMAL  FREQUENCIES 

.00  a CC  0.0 

• 014 

3000  .2326000  .2443187 

_ .2716625 

.3146312  .3615062  .4161937  .4708812 


Figure  27  - Frequency  Response  for  10  Stage  1:2  Digital  Filter 


Frequency  Response  for  55  Stage  1:1  Digital  Filter 
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Figure  32  - Frequency  Response  for  16  Stage  1:4  Digital  Filter 


Frequency  Response  for  10  Stage  1:2 
Digital  Filter  Rounded  off  to  12  Bits 


Figure  34  - Frequency  Response  for  10  Stage  1:4 

Digital  Filter  Rounded  off  to  12  Bits 


Frequency  Response  for  25  Stage  1:4 
Digital  Filter  Rounded  off  to  12  Bits 
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4.2  Hacro  Instructions 

A nacro-lastructlon  Is  a set  of  ones  and  zeroes  that  instructs  the  machine 
on  processing  tasks  to  be  performed.  From  the  algorithm.  It  can  be  seen  that 
the  filtering  process  is  the  most  Important  task  that  the  NBSCU  has  to  perform. 
There  are  several  pieces  of  Information  that  have  to  be  supplied  to  the  machine 
for  the  filtering  operation.  These  Include  the  following: 

a.  Starting  Address  - The  machine  has  to  be  instructed  as  to  where  to 

fetch  the  coefficients  for  convolution. 

b.  Ending  Address  (or  number  of  coefficients)  - The  machine  has  to  know 

when  the  convolution  ends  so  it  can  stop.  The 
two  specifications  (ending  address,  or  number  of 
coefficients)  are  equivalent.  However,  due  to 
ease  of  implementation,  one  specification  may  be 
preferred  to  the  other. 

c.  Resampling  Rate  - The  machine  has  to  know  how  often  the  samples  are 

to  be  processed.  For  a 4:1  filter,  since  only 
every  fourth  sample  Is  needed,  the  other  three 
output  samples  may  not  be  computed.  This  will 
represent  a significant  reduction  in  loading. 

Based  on  the  requirements  on  the  filtering  operation,  a macro-instruction 
format  Is  obtained  and  is  shown  in  figure  39. 

It  Is  broken  Into  6 sub-fields  with  the  following  bit  assignments: 

Bit  Assigned  Function 

0 This  is  the  bit  that  determines  whether  a single  channel  (0) 
algorithm  (as  in  AM  transmit)  or  a double  channel  (1)  algorltlm 
(as  in  AM  receive)  Is  used. 

1 This  bit  determines  the  rate  at  which  the  memory  counters  count. 

For  a dual  channel  operation,  the  rate  Is  half  of  that  of  a 
single  channel  operation  because  two  multiplications  are  required. 

2 This  bit  is  reserved  for  single  sideband  transmit  mode. 

3 This  bit  represents  a frequency  translation  process  occurring 
after  a filtering  process. 

4-13  These  ten  bits  represent  the  starting  address  for  the  filter 
coefficients. 

14-23  These  ten  bits  specify  the  final  location  for  the  filter 
coefficient. 
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Figure  39  - Macrolnetructlon 
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24-31  These  are  the  bits  chat  specify  the  resampling  rate  and  are  defined 
as  follows: 

(a)  If  Che  sample  rate  remains  unchanged  a filter  resampling 
specification  of  11  Is  used. 

(b)  If  Che  sample  rate  Increases  after  filtering  (e.g.,  four  times 
Che  rate)  a specification  of  14  Is  used. 

(c)  If  Che  sample  rate  decreases  after  filtering  (e.g.,  one  fourth 
Che  rate)  a specification  of  41  Is  used. 

In  all  Che  above  cases,  the  first  digit  Is  represented  by  the  first 
four  bits  In  the  field  while  Che  second  digit  Is  represented  by  the 
last  four  bits  In  the  field. 

32-35  These  four  bits  represent  Che  function  to  be  performed  after  Che 
present  stage  of  filtering  and  appropriately  conditions  the  bits 
for  use  In  Che  next  stage  of  operation.  These  bits  are  defined  as 
follows: 

0000  This  code  calls  for  Che  arctangent  demodulation  routine  In 
Che  AM  case. 

0001  This  code  calls  for  Che  arctangent  demodulation  routine  In 
Che  FM  case. 

0101  This  Is  the  frequency  translation  routine. 

0100  This  Is  used  for  pre-emphasis  filtering  and  frequency 

translation  that  may  be  Incorporated  Into  the  algorithm  later 
on. 

1000  This  Is  for  de-emphasls  filtering  and  will  be  using  In  con- 
junction with  Che  pre-emphasis  filter. 

1100  This  Is  for  AM  transmit  In  which  a carrier  Is  added  to  the 
waveform  as  well  as  frequency  translation. 

1111  This  Is  used  for  consecutive  filtering  process. 

Therefore,  If  a 16  stage  filter  Is  to  be  used  In  a single  channel  mode 
with  Che  starting  address  In  Che  coefficient  listing  located  at  address  10, 
and  assuming  a 1:1  resampling  rate,  Che  following  macrocode  can  be  obtained 
for  this  process. 

001100000010010000010101000100011111 


4.3  The  Arctangent  Demodulator 


The  Arctangent  algorithm  Is  a fast  and  efficient  algorithm  to  determine 
the  phase  and  amplitude  of  a complex  number.  The  process  Is  essentially  a 
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rectangular  to  polar  coordinate  conversion  process.  Given  a complex  number, 
which  can  be  represented  by  the  vector  diagram  In  figure  40,  the  amplitude  can 
be  obtained  by  using  the  equation 

A ■ R cos  0 + Q sin  0 . 

The  angle  0 Is  given  by: 

0 ■ tan"^  ^ 

R 

Notice  that  regardless  of  the  angle  0 , the  amplitude  Is  always  positive. 

This  Is  equivalent  to  the  diode  detector  output  of  an  AM  demodulator  which  Is 
always  positive.  Only  two  multiplications  and  one  addition  are  required  to 
obtain  the  amplitude  provided  that  the  phase  Is  known.  It  Is,  therefore, 
obvious  that  as  far  as  the  demodulation  process  goes,  phase  demodulation  Is 
the  most  efficient.  If  the  Arctangent  Demodulator  Is  to  be  used. 

Obtaining  the  absolute  phase  angle  of  the  vector  Is  a little  difficult 
because  It  Involves  a division  process  as  well  as  an  arctangent  process. 
Straight  Implementation  of  the  division  process  and  the  series  method  of 
obtaining  the  arctangent  of  a number  Is  very  Inefficient  and  involves  many 
multiplications.  These  two  methods  are  therefore  abandoned.  A more  efficient 
way  to  Implement  the  two  functions  Is  by  means  of  table  look  up.  By  putting 
two  tables  (I/X,  arctangent)  on  read  only  memories  (ROMS) , the  phase  angle 
can  be  obtained  by  merely  one  multiplication.  Therefore,  four  tables,  as  seen 
from  the  equation,  are  required  for  the  Arctangent  Demodulator  for  both  ampli- 
tude and  phase  computation.  These  are  the  1/X,  tan~l,  sine  and  cosine  tables. 
Since  the  sine  and  cosine  can  be  obtained  from  one  another,  one  table  can  be 
eliminated.  Furthermore,  since  the  value  of  the  sine  and  cosine  function 
are  unique  up  to  90®,  the  required  sine  table  will  be  one  fourth  of  the  full 
sine  table.  To  reduce  the  size  of  the  arctangent  table,  the  magnitude  of  the 
two  vector  components  can  be  obtained  and  compared.  The  larger  of  the  two 
vectors  can  be  used  to  address  the  1/X  table.  Multiplication  of  the  table 
output  and  the  smaller  vector  will  produce  a number  that  Is  less  than  1,  cor- 
responding to  45*.  Therefore,  a tan“^  table  ranging  from  0 to  45®  Is  all  that 
Is  required.  The  procedure  for  obtaining  the  phase,  and  amplitude  Is  as 
follows : 

1.  Obtain  magnitude  of  vectors  and  store  original  signs  of  vectors. 

2.  (^pare  the  two  numbers  and  use  the  larger  of  the  two  to  address 
the  1/X  cable. 

3.  Multiply  Che  output  of  Che  1/X  cable  with  the  smaller  number  to  form 
Che  product  |y|/]x|. 

4.  Use  Che  product  to  address  Che  arctangent  ROM  to  obtain  the  angle 
In  Che  first  octant. 

For  getting  Che  magnitude  of  the  vector,  Che  angle  can  be  put  Into  Che 
proper  quadrant  and  the  corresponding  values  of  Che  sine  and  cosine  of  the 
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angle  In  the  first  quadrant  are  obtained  from  the  sine  table.  The  amplitude 
Is  equal  to  R cos  0 + Q sin  0.  Since  the  amplitude  Is  always  positive,  and 
also  cos  0 and  sin  0 Is  positive  In  the  first  quadrant,  the  following  equation 
can  be  used  without  affecting  the  result. 

/R/  cos  0 + /Q/  sin  0 

where  0°  1 0 1 90°  and  Is  the  angle  that  has  the  same  absolute  sine  and  cosine 
value  as  the  original  angle.  The  absolute  phase  angle  can  then  be  determined 
by  using  table  12.  Notice  that  the  angle  ranges  from  +180°  to  -180°  because 
bipolar  2's  complement  arithmetic  is  used. 


TABLE  12  - PHASE  DETERMINATION  OF  ARCTANGENT  DEMODULATOR 


_R_ 

1r1  - loi  Octant 

Absolute  Angle 

+ 

+ 

+ 

1 

Tan“l 

|Q|/|R| 

+ 

+ 

2 

1 

O 

O 

Tan-1  [R1/!q| 

- 

+ 

3 

90°  + 

Tan-1  IrI/|q| 

- 

+ 

+ 

4 

180°  - 

Tan-1  |q|/|r1 

- 

- 

+ 

5 

-180°  + 

Tan-1  |Q|/|r| 

- 

- 

i 

-90°  - 

Tan-1  1r1/1q1 

+ 

- 

7 

-90°  + 

Tan-1  |R|/|q| 

+ 

- 

+ 

8 

-Tan-1 

!q1/|r| 

NBSCU  Design  Requirements 

Two  of  the  most  important  functions  that  the  narrowband  signal  conversion 
unit  has  to  perform  Is  that  of  filtering  and  demodulation.  Since  the  unit  has 
to  handle  FM  signals  which  are  very  phase  sensitive,  non-recurslve  finite  im- 
pulse response  (FIR)  filters  were  chosen  over  Infinite  impulse  response  (HR) 
filters.  In  the  FM  and  SSB  mode,  notice  that  the  output  of  one  filter  Is  the 
Input  to  another  filter,  and  therefore,  the  filtering  circuits  must  have  the 
capability  to  loop  the  output  back  to  the  Input  of  the  filter  to  be  reprocessed. 
Furthermore,  since  the  filters  are  performing  the  resampling  and  desampling 
functions,  this  Imposes  a condition  on  the  Input  and  output  buffer  size.  The 
maximum  requirement  occurs  In  the  single  sideband  receive  and  transmit  mode 
In  which  32  Input  samples  are  required  to  produce  one  output  sample.  Thus, 
the  buffers  used  must  be  at  least  32  words  In  length.  A clock  generator  Is 
also  required  to  supply  the  required  clocks  to  the  A/D  and  D/A  converters. 

A multiplier,  which  Is  the  most  Important  element  In  digital  filtering,  must 
be  capable  of  operating  fast  enough  to  do  the  required  number  of  multiplications 
In  real  time.  The  various  fimctlons  that  are  to  be  performed  should  be  micro- 
processor controlled.  In  the  event  that  the  microprocessor  falls.  It  Is 
desirable  that  a standby  function  can  be  loaded  manually  by  the  operator. 

This  Is  also  one  of  the  redundant  features  In  the  TIES  program. 
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5.  NBSCU  CIRCUIT  DESIGN 


Circuits  were  designed  based  on  the  functions  that  the  narrowband  signal 
conversion  unit  has  to  perform  and  can  be  separated  Into  two  categories  — 
analog  and  digital.  A top  level  block  diagram  for  the  narrowband  signal  conver- 
sion unit  Is  shown  In  figure  41. 

Discussions  on  the  analog  and  digital  circuits  are  given  in  the  following 
sections. 

5.1  Analog  Interface  Circuits 

The  analog  circuits  serve  as  the  Interface  between  the  analog  world  and 
the  signal  processor.  The  Interface  consists  of  a somple-and-hold  circuit,  an 
A/D  converter  at  the  Input  and  D/A  converter  at  the  output.  A circuit  that 
generates  the  sampling  rate  suitable  for  conversion  Is  also  part  of  the  analog 
Interface  circuit  which  Is  shown  In  figure  42.  All  circuit  components  are 
commercially  available  hardware. 

Upon  command  from  the  sampling  rate  generator,  the  Input  signal  to  the 
sample-and-hold  circuit  Is  sampled.  The  sampled  value  Is  then  held  by  a capl- 
citor.  The  A/D  convertor  will  produce  a 10  bit  word  In  two's  complement  nota- 
tion. The  sampling  rate  generator  Is  simply  a programmable  dlvlde-by-N  circuit 
that  takes  In  a frequency  of  230.4  KHz  and  produce  a sub-multiple  of  this 
frequency.  Depending  on  whether  the  transmit  or  receive  mode  is  used,  the 
higher  or  lower  frequency  Is  routed  appropriately  to  the  A/D  or  D/A  converter. 
Notice  that  no  phase  lock  loop  Is  required  to  generate  the  sampling  frequencies 
because  the  resampling  and  desampling  filters  are  multiples  or  reciprocals 
of  the  Input  sampling  rate. 


5.2  Digital  Circuits 

A detailed  block  diagram  for  the  NBSCU  Is  shown  In  figure  43. 

The  digital  circuit  can  be  partitioned  into  two  parts  — the  loading 
circuits  used  to  load  the  filter  coefficients  from  mass  storage  memory  Into  the 
memory  of  the  signal  processor  and  the  circuits  that  process  the  input  samples. 

5.2.1  Loading  Circuits 

The  loading  circuit  la  shown  In  figure  44.  The  operation  of  this  circuit 
Is  described  as  follows: 

Toggle  switch,  SOI,  In  the  open  position,  enables  EPROMS  1 and  2 and  the 
B Inputs  to  MUl.  By  depressing  the  preset  switch  S02,  the  appropriate  flip- 
flops  and  counters  are  preset.  Depressing  switch  S03  causes  DOlOl  to  change 
state,  thereby  turning  the  master  clock  on  at  A020103.  The  clock  then  goes 
through  a counter  which  has  outputs  that  are  sub-multiples  of  the  basic  clock 
rate.  This  Is  necessary  because  the  EPROMS  have  a cycle  time  that  Is  much 
longer  than  the  200  nsec  master  clock  period.  Using  the  divide  by  lb  output 
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as  the  clock  to  the  EPROM  address  counter,  the  contents  of  the  first  location 
(which  consists  of  Che  number  of  stages  of  filtering  involved  in  the  Dod/deno 
process)  is  loaded  into  flip-flops  D04,  DOS  while  the  content  of  location  2 
(which  contains  the  required  number  of  constants  and  coefficients)  is  loaded  into 
flip-flops  D06,  D07.  The  36  bit  macro  instructions,  each  of  which  occupy  5 loca- 
tions in  Che  EPROM  are  loaded  into  the  instruction  RAMs  which  are  addressed  by 
counters  BOS,  B06.  Comparators  C04,  COS,  which  compare  Che  number  of  stages 
required  for  the  mod/demod  process  and  the  counter  outputs,  produce  a signal 
when  all  the  macro  Instructions  are  loaded  into  the  instruction  RAMs.  The 
coefficient  RAMs  then  begin  to  fill  up.  The  addressing  device  in  this  case 
is  a set  of  arithmetic  logic  units  because  of  Che  resampling  requirement  idilch 
will  be  seen  later  on.  Upon  completion,  comparators  COl,  C02,  C03  produce  a 
signal  Chat  changes  Che  state  of  flip-flop  D0201,  which  turns  a LED  on  indica- 
ting Chat  Che  processor  is  loaded. 


5.2.2  Signal  Processor  Circuits 


Upon  completion  of  Che  loading  process,  the  signal  processor  is  left 
in  Che  idle  state.  It  will  reoialn  in  Che  idle  state  until  enough  samples  re- 
quired for  processing  are  gathered.  The  number  of  samples  required  can  be 
determined  by  Che  overall  resampling  rate.  The  circuit  for  buffer  loading  is 
shown  in  figure  45. 

When  Che  required  samples  are  acquired,  a signal  is  produced  at  0240308 
and  shaped  by  the  one  shot  circuit  MX02  which  produces  a 30  nsec  pulse  that 
changes  the  state  of  D0802  which  turns  00401  on.  The  signal  processor  now 
becomes  active  with  Che  master  clock  turned  on  at  0040103.  Three  reset  pulses 
are  produced  and  designated  RESET  1,  RESET  2,  RESET  3.  RESET  1 is  used  Co  re- 
set devices  that  need  Co  be  reset  after  all  the  input  samples  are  processed. 
RESET  2 is  used  Co  reset  devices  Chat  need  to  be  reset  after  each  filtering 
process.  RESET  3 is  used  to  reset  devices  that  need  to  be  reset  after  each 
convolution. 

After  Che  proper  devices  are  reset  by  the  three  reset  pulses,  Che  input 
data,  which  are  being  held  by  Che  input  buffer  begins  the  transfer  Co  the  high 
speed  RAMs.  The  data  transfer  is  complicated  due  to  the  resampling  process. 

The  circuit  is  shown  in  figure  46. 

Inputs  to  comparators  COl  and  C02  will  be  coming  from  the  macro  Instruc- 
tion RAMs  representing  Che  resampling  rate  used  for  Che  particular  filter. 

Since  a resampling  filter  of  4:1  requires  the  processing  of  1 sample  out  of 
4 input  samples,  four  samples  are  loaded  into  the  high  speed  RAMs.  However, 

Che  convolution  process  only  cakes  place  after  Che  fourth  sample  is  loaded. 

On  the  other  hand,  if  a 1:4  filtering  function  is  to  be  performed  all  input 
samples  are  used  (including  the  zeroes  used  for  zero  filling).  Therefore, 

Che  inputs  to  comparator  COl  determine  the  overall  samples  required  to 
produce  one  sample.  However,  Che  input  Co  Che  comparator  C02  does  not 
represent  the  number  of  zeroes  to  be  filled.  Instead,  the  nximber  of  zeroes 
needed  is  one  less  than  the  input  to  comparator  C02.  Provision  is  made  so 
Chat  Che  state  of  the  counter  BOl  can  be  stored  in  memory  R05  in  case  not 
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enough  samples  are  available  for  processing.  When  enough  samples  are  loaded 
in  the  high  speed  RAMs  for  one  convolution,  a signal  from  comparator  COl  will 
turn  on  the  clock  to  the  convolutional  filter  at  0220102. 

Notice  that  enough  samples  must  be  loaded  Into  the  high  speed  RAMs  for 
convolution  to  begin.  The  processor  will  go  to  the  next  macro  Instruction  If 
the  samples  are  Insufficient,  l.e.,  the  program  counters  BOS,  B06  will  increment 
when  this  condition  occurs  even  no  convolution  takes  place.  This  operation  Is 
controlled  by  the  circuit  shown  on  figure  47. 

The  high  speed  RAMs  are  addressed  by  counters  B09,  BIO,  Bll,  which  are 
loaded  with  an  Initial  address  for  the  particular  filter  and  Increment  each  time 
new  data  Is  deposited  onto  the  RAMs.  Decrementation  begins  when  actual  con- 
volution takes  plac.t.  The  address  at  which  the  counters  stop  counting  Is 
stored  In  RAMs  ROl,  R02,  R03  because  these  data  have  to  be  retrieved  for 
processing  the  next  set  of  data.  This  Is  shown  In  figure  48. 

Since  each  time  a new  sample  is  loaded  Into  the  high  speed  RAM,  the 
address  counters  Increment  by  one  to  get  ready  for  the  next  Input,  these 
co«inters  have  to  be  decremented  by  one  when  convolution  occurs.  For  a parti- 
cular filter  operation,  the  beginning  and  ending  addresses  for  the  coefficients 
are  used  and  the  circuit  must  have  the  capability  to  reset  the  counter  to  the 
starting  address  during  the  loading  process  and  must  be  able  to  reset  to  the 
ending  address  during  the  convolution  process.  This  Is  accomplished  by  en- 
abling the  A Input  of  multiplexers  M29,  M30,  M31  to  detect  the  condition  that 
the  ending  address  is  reached.  The  M24,  M2S,  M26,  M27,  M28  multiplexers  are 
then  enabled  so  that  the  beginning  address  lines  are  used  for  loading  Into 
counters  B09,  BIO,  Bll.  At  the  beginning  of  the  filtering  ftinctlon,  the 
previous  address  Is  loaded  Into  the  counters  by  RESET  3 from  RAMs  ROl,  R02, 

R03.  Convolution  begins  upon  completion  of  the  loading  process. 

The  circuit  for  the  convolution  process  Is  shown  In  figure  49. 

Two  shift  registers  SOI,  S02  are  used  for  controlling  the  timing  re- 
quirements for  the  two  multipliers,  the  ALUs  and  the  ALU  output  registers. 

SOI  Is  for  the  real  channel  while  S02  Is  for  the  quad  channel  in  the  two  channel 
operating  mode.  Data  from  the  two  channels  are  alternately  loaded  Into  the 
Input  registers  of  the  multiplier  while  outputs  of  the  ALUs  are  loaded  appropri- 
ately Into  the  corresponding  accumulator  registers.  In  the  single  channel 
operating  mode,  the  multiplier  circuitry  Is  the  same  as  In  the  dual  channel 
m^e.  However,  data  now  comes  from  one  channel  at  all  times  and  only  one 
accumulator  register  la  activated.  To  operate  the  circuit  at  the  most 
efficient  rate,  the  address  counters  B09,  BIO,  Bll  have  to  be  clocked  at 
twice  the  rate  so  that  the  two  multipliers  can  operate  at  the  same  rate  as  In 
the  two  channel  case.  This  Is  accomplished  by  flip-flops  D2S02  upon  command 
from  the  Instruction  RAM.  Notice  that  In  both  the  single  and  dual  channel 
cases  after  the  last  data  word  Is  loaded  Into  the  multiplier,  three  clock 
pulses  are  required  before  data  Is  available  from  the  multiplier.  This  Is 
taken  care  of  by  counter  B20.  Signals  coming  out  of  counter  B20  will  trigger 
the  state  of  flip-flop  DllOl  to  change  Indicating  that  convolution  Is  completed. 
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Addressing  Circuit  for  High  Speed  RAMs 


Figure  49  - Control  Circuit  for  Convolutional  Filtering 


Figure  40  - Real  and  Quad  Components  of  a Complex  Number 
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The  control  circuits  for  the  high  speed  RAMs,  the  multipliers,  the  ALUs  and 
associated  registers  are  shown  In  figures  SO  through  55. 

Depending  on  the  process  to  be  performed  after  filtering,  the  data 
stored  in  the  accumulator  registers  have  to  be  conditioned.  This  can  be  done 
quite  easily  by  a set  of  shift  registers.  This  Is  shown  In  figure  56. 

The  change  In  state  of  flip-flop  DllOl  enables  the  1 of  16  decoder  to 
which  the  last  four  bits  of  the  macro-instruction  are  wired.  The  1 of  16  de- 
coder, when  enabled,  activates  the  proper  set  of  shift  registers.  As  the  zero 
logic  level  propagates  down  the  shift  register  the  corresponding  functions  are 
accomplished.  During  these  processes,  the  operation  of  the  multipliers,  the 
ALUs,  as  well  as  the  accumulator  registers  may  have  to  be  activated  again. 

The  control  signals  for  these  functions  can  be  found  In  figures  50,  52,  54, 
respectively. 

The  various  tasks  to  be  performed  after  filtering  are  as  follows: 

a.  Filtering 

This  Is  shown  In  figure  57  and  gives  the  capability  to  the  proces- 
sor to  perform  cascade  filter  operation.  The  data  are  not  changed  and  are  loaded 
back  into  the  real  and  quad  RAMs  appropriately.  To  simplify  circuit  design, 
both  real  and  quad  accumulators  are  used  even  In  the  single  channel  mode  In 
which  only  one  register  Is  used. 

b.  AM  Transmit  - Modulation 

After  the  last  stage  of  filtering  In  the  AM  transmit  mode,  a carrier 
Is  added  to  the  filter  output  sample  and  the  resulting  signal  Is  mixed  up  to 
57.6  KHz  for  transmission.  This  circuit  has  to  add  to  the  accumulator  register 
a constant  representing  the  carrier  and  multiplying  the  resulting  output  by 
1,0, -1,0 — the  four  samples  of  a sinusoidal  wave  for  frequency  trans- 

lation. Since  the  ALUs  can  perform  a negation  function,  no  multiplication  Is 
actually  carried  out.  Furthermore,  since  multiplication  by  zero  produces  a 
zero  output,  only  every  other  output  sample  Is  required  to  be  computed.  There- 
fore, for  this  particular  process,  the  filter  preceding  the  carrier  addition 
process  can  be  set  to  be  a 2:4  filter.  This  will  reduce  the  overall  number  of 
multiplication  required  for  the  process.  Notice  that  the  specification  Is 
2:4  and  not  1:2.  A 2:4  specification  Indicates  that  2 samples  are  required 
to  produce  1 output  sample  and  that  three  out  of  four  input  samples  are 
zeroes,  which  is  what  the  normal  input  would  be  If  a 1:4  resampling  were  to 
be  done.  A 1:2  specification  will  produce  one  output  sample  for  every  Input 
sample  which  has  zero  value  at  every  other  sample.  This  Is  Incorporated 
Into  the  circuit  shown  In  figure  57. 

c.  AM  Receive  - Demodulation 

The  circuit  diagram  for  this  case  Is  shown  In  figure  58. 


i.uiiccion  can  be  loaded  manually  by  the 

This  is  also  one  of  the  redundant  features  in  the  TIES  program. 


operator. 
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Figure  50  - High  Speed  RAMs  Control  Circuit 


Hode  Control  Circuit  for  Arlthmatic  Logic  Units 


Control  Circuit  for  Arlthmatlc  Logic  Units 
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Reprocessing  of  Filtered  Data  for  AM  Transmit,  SSB  Transmit /Receive 
and  Arctangent  Demodulation  for  FM  Receive  Modes 
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The  function  of  this  circuit  Is  to  obtain  the  amplitude  of  a vector 
given  In  rectangular  coordinates  and  In  the  arctangent  process  described 
earlier.  Notice  that  the  absolute  angle  Is  not  computed  since  the  amplitude 
Is  always  positive. 

d.  FM  Transmit  - Modulation 

The  circuit  for  FM  modulation  Is  shown  In  figure  59. 

The  output  of  the  1:4  filter  Is  integrated  to  give  a phase  represen- 
tation of  the  Incoming  signal  which  Is  used  to  address  a slne/coslne  table  to 
produce  the  sine  and  cosine  value  of  the  angle.  The  frequency  translation 
process  Is  obtained  by  changing  the  signs  of  these  numbers  according  to  the 
FM  transmit  algorithm. 

e.  FM  Receive  - Demodulation 

The  circuit  for  the  FM  receive  demodulation  process  is  Incorporated 
In  figure  58.  This  Is  the  same  arctangent  process  used  In  the  AM  case.  In 
the  case  for  FM,  the  absolute  phase  Is  computed  and  Is  to  be  differentiated 
to  obtain  the  frequency  of  the  signal. 

f.  SSB  Transmit  and  Receive 

FollotdLng  the  last  stage  of  filtering  in  the  single  sideband 
receive  or  transmit  mode,  a frequency  translation  can  be  seen  from  the  algo- 
rithm. This  la  similar  to  the  FM  transmit  case  In  which  only  every  other 
sample  from  the  two  channels  are  needed  to  produce  the  result.  This  Is  in- 
corporated Into  figure  57.  The  frequency  translation  In  this  case  Is  equivalent 
to  multiplying  the  data  by  +1,  +1,  -1,  -1 etc. 

Notice  that  the  coefficient  and  constant  RAMs  which  are  normally 
addressed  by  the  ALUs  have  to  be  addressed  differently  depending  on  which 
constant  Is  required  during  the  demodulation  process.  This  circuit  diagram 
Is  shown  In  figure  60  which  Is  also  the  control  logic  for  addressing  the  RAM 
In  the  loading  diagram. 

After  processing  one  set  of  samples,  the  resulting  output  can  be  D/A 
converted  back  Into  an  analog  signal  or  the  data  can  be  put  back  into  menory 
for  further  filtering.  For  reasons  of  simplicity  In  control  circuitry,  two 
sets  of  RAMs  are  needed.  While  the  contents  of  one  set  of  RAMs  are  used  as 
Inputs  to  the  processor,  the  other  set  stores  the  processed  output.  The 
roles  are  reversed  between  the  two  sets  of  RAMs  after  processing  one  set  of 
data.  To  accommodate  the  dual  channel  case,  therefore,  four  sets  of  RAMs 
are  used.  However,  the  address  counters  remain  as  two  sets.  To  eliminate 
the  time  required  for  transferring  the  Input  samples  to  these  RAMs,  the  input 
buffer  Is  chosen  to  be  the  other  set  of  the  same  RAM.  By  ping  ponging  this 
set  of  RAMs  with  one  of  the  two  sets  of  RAMs  In  the  real  channel,  data  trans- 
fer can  be  accomplished  readily.  The  control  circuits  can  be  seen  on  figure 
45  and  the  Input/output  wiring  can  be  seen  on  figure  61. 
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HARDWARE  REQUIREMENTS 

To  Implement  the  circuits  in  the  design,  high  speed  ICs  have  to  be  used 
because  of  the  real  time  requirement.  Among  the  different  processes  Involved, 
multiplication  is  most  critical  and  time  consuming.  The  two  candidates  that 
give  adequate  multiplication  performance  will  be  the  AMD  25S0S  multiplier  and 
also  the  TRW  SJ  series  multiplier.  The  AMD  device  has  a very  slight  edge 
over  the  TRW  multiplier  In  speed.  However,  the  AMD  device  only  comes  as  a 
2X4  bit  package  and  would  Involve  eighteen  24-pln  ICs  plus  5 high  speed 
carry  look  ahead  ICs  to  realize  a 12  bit  by  12  bit  multiplier.  Because  of 
the  real  estate  Involved,  the  AMD  device  was  abandoned.  The  TRW  device  Is 
an  LSI  device  with  64  pins  and  can  perform  multiplication  at  a rate  of  175 
nsec,  and  draws  a maxlnnun  of  850  ma  equivalent  to  4.25  watts  at  5 volts. 

All  ALUs  are  Schottky  ICs.  For  a 24  bit  addition,  with  high  speed 
carry  look  ahead  circuits,  the  addition  time  Is  only  19  nsec  and  presents  no 
problem  of  loading  the  circuit  down. 

The  overall  constraint  on  slowing  the  operating  speed  comes  from  the 
address  counters  B09,  BIO,  Bll  as  well  as  the  RAMs.  This  Is  especially  true 
In  the  single  channel  mode  In  vdilch  the  counters  are  to  be  Incremented  at 
twice  the  dtial  channel  mode.  The  worst  propagation  delay  for  a programmable 
up/down  counter  Is  36  nsec,  vdiile  the  address  access  time  for  the  RAMs  Is 
30  nsec,  making  a total  of  138  nsec.  This  Is  equivalent  to  a clock  rate  of 
about  7 MHz.  To  operate  at  this  rate,  two  multipliers  are  required.  This 
will  give  a safety  margin  as  well  as  the  capability  for  expansion  in  the 
future. 


6.  NBSCU  Specifications 

The  following  are  specifications  for  major  components  needed  as  well  as 
various  design  parameters  based  on  the  present  NBSCU  design. 

6.1  SAMPLE  AND  HOLD 

Manufacturer:  DATEL 

Model  No.:  SHM-UH 

Input  Range:  ±5V  FS 

Input  Impedance:  100  Megohms  shunted  with  20  pf. 

Sample  command:  35  t 10  nsec 

0 to  5 volts  positive  pulse 
3 nsec  max.  rise  and  fall  time 

Sample  command  Input  Impedance:  50  ohms  - requires  100  ma  of  source 

current 
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Output  Voltage  Range:  15V  FS 
Output  Current:  130  ma  max. 

Output  Impedance:  10  ohms,  max. 

Acquisition  Time:  35  nsec  for  full  scale  Input 
Aperture  Time:  200  psec. 

Output  Slewing  Rate:  500V/ysec 
Maximum  Sample  Rate:  10  MHz 
Hold  Decay  Rate:  50  uV/ysec 

6.2  A/D  Converter 
Manufacturer : DATEL 
Model  Number : ADC-G3C 

Output  Coding:  2's  complement  (reverse  coding) 

No.  of  Bits:  10 

Analog  Input  Range:  ±5  volts 

Input  Impedance : 2IS2 

Start  Conversion:  2V  min.  to  7V  max.  positive  pulse. 

30  nsec  min.  width 

"1"  - Resets 

"0"  “ Starts  Conversion 

End  of  Conversion  Output:  0 ■ Conversion  complete 

1 ■ Reset  and  Conversion  Period 

Conversion  Rate:  1 MHz 

6.3  D/A  Converter 
Manufacturer : DATEL 
Model  No.:  DAC-VR  12B3C 
Input  Coding:  2's  complement 
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Input  Loading:  1 TTL  load 

Data  Strobe:  Information  must  be  present  at  the  register  Inputs  of 

the  D/A  converter  prior  to  strobing.  Transfer  of  Informa- 
tion to  the  output  pins  occurs  when  the  strobe  Input  goes 
from  1 to  0. 

Update  Rate:  10  MHz  max. 

Output  Voltage:  ±5  volts 

Output  Settling  Time:  <2  usee  to  ± 0.025Z  of  FS 
6.4  Processor 

Input  Buffer  Size:  256,  12  bit  words 

Scratch  pad  memory  size:  1,024,  12  bit  words 

Macroinstruction  memory  size:  Sixteen  36  bit  words 

Multiplication  Rate:  10  MHz 

Addition  Rate:  SO  MHz 

Maximum  Sampling  Rate:  5 KHz  - 500  KHz 

No.  of  IC's  used:  360 

Power  Supply:  5 volt,  -12  volts 

Current  drawn:  20  amperes  at  5 volt 
70  ma  at  -12  volts 

Output  Buffer  Size:  64,  12  bit  words  (FIFOs) 

No.  of  Circuit  Boards:  One  19"  x 8"  Multiplier  - ALU  Board 

One  19"  X 8"  Memory  Board 
One  10"  X 10"  I/O  and  Control  Board 
One  8"  X 8"  Analog-Digital  Conversion  Board 
These  boards  are  shown  In  figure  62. 
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APPENDIX  A 


NOISE  SOURCES 


In  Che  course  of  digital  processing  of  signals,  there  are  many  places 
Chat  noise  may  be  introduced  chat  may  affect  the  qu^lty  of  Che  signal.  In 
the  following  paragraphs,  some  of  these  noise  sources  are  investigated. 

A.l  Noise  Sources  from  a Sampling  Process 

The  quality  of  the  signal  begins  Co  deteriorate  the  moment  it  is 
sampled  — assuming  sampling  is  high  enough  in  frequency  so  Chat  there  is  no 
aliasing  effect.  The  first  group  of  noise  sources  come  from  the  sample-and- 
hold  circuit  which  freezes  Che  fast  moving  signal  for  the  A/D  converter  to 
make  the  conversion.  Those  noise  sources  are  listed  below: 

1)  OFFSET  - This  is  the  extent  Co  which  Che  output  deviates  from  zero 

with  zero  input  and  is  usually  a function  of  time  and  tem- 
perature and  can  be  adjusted. 

2)  NONLINEARITY  - This  is  the  amount  by  which  the  output  versus  input 

deviates  from  a straight  line. 

3)  SCALE  FACTOR  ERROR  - This  is  Che  aiaount  by  which  the  output  deviates 

from  specified  gain  (usually  unity) . This  para- 
meter can  also  be  adjusted. 

4)  APERTURE  UNCERTAINTY  - This  is  the  time  elapsed  between  the  command 

Co  hold  and  the  acttial  opening  of  Che  hold 
switch. 

5)  DROOP  - This  is  the  drift  of  Che  output  caused  by  the  discharge  of 

the  storage  capacitor. 

Among  Che  above  noise  sources  from  the  sample-and-hold , the  most  impor- 
tant ones  are  the  aperture  uncertainty  and  Che  droop. 

The  aperture  uncertainty  has  two  components  — a nominal  time  delay,  and 
an  uncertainty  caused  by  Jitter  or  variation  from  time  Co  time  and  unit  to 
unit  and  is  shown  in  figure  Al. 

The  aperture  uncertainty  of  a sample-and-hold  unit  dictates  the  rate  that 
Che  input  signal  can  vary  given  chat  the  signal  is  to  be  resolved  to  a certain 
desired  minimum  value  as  in  Che  case  of  an  A/D  converter.  Assuming  that  a 
signal  of  the  form  A sin  cot  is  sampled  by  a sample-and-hold  circuit.  The 
rate  at  which  the  signal  varies  with  time  will  be  Che  time  derivative  of  the 
signal  and  is  given  by  the  following  equation: 

^A  sin  cot  ■ (oA  cos  cot  volCs/sec 


A-1 


NADC-79153-40 


If  an  A/D  converter  is  to  be  used  to  put  the  sampled  value  Into  10  binary 
bits,  the  aperture  error  should  not  cause  the  output  to  vary  more  than  1/2  LSB. 
Assuming  the  A/D  converter  has  1 volt  as  Its  full  scale,  therefore  1/2  LSB  ■ 2~il 
volts.  The  maximum  rate  at  which  the  input  signal  can  vary  is  (uA  volts/sec, 
and  therefore,  the  maxlaami  aperture  error  is  2~H/a)A  seconds.  Therefore,  to 
specify  the  aperture  error,  one  must  know  the  weight  of  the  LSB  of  the  A/D 
converter  and  the  strength  of  the  incoming  signal.  Notice  that  if  oiA  becomes 
larger  indicating  that  the  frequency  is  larger  (A  is  a constant) , a smaller 
aperture  uncertainty  is  required. 

As  aperture  time  is  important  to  fast  moving  signals,  the  droop  is 
Important  to  slow  moving  waveforms.  A typical  diagram  for  the  droop  is  shown 
in  figure  A2  and  is  the  effect  caused  by  the  discharge  of  the  capacitor  in  the 
hold  circuit.  Since  slow  moving  signals  Imply  a slower  sampling  rate,  the 
capacitor  will  have  more  time  to  discharge  than  in  higher  rate  sampling  cir- 
cuits. Therefore,  given  the  LSB  for  the  A/D  converter  using  in  conjunction 
with  the  sample-and-hold  with  a certain  droop,  the  minimum  sampling  rate  can 
be  determined.  For  example,  assume  a sample  and  hold  has  a droop  rate  of 
50  yv/usec  and  a 1 volt  full  scale  10  bit  A/D  converter  has  2~H  volt  as  the 
weight  of  half  its  LSB.  The  minlaium  sampling  rate  should  be: 

f ■ [<2~H  volt)(l  sec/50  volt)]"^  ■ 102.4  KHz 

Sampling  less  than  102.4  KHz  may  yield  an  erroneous  result  at  the  output 
of  the  A/D  converter. 

The  source  of  error  coming  from  the  A/D  converter  is  known  as  quantization 
noise.  The  noise  arises  because  of  the  fact  that  only  a finite  number  of  bits 
are  used  to  represent  an  analog  voltage.  This  is  known  as  truncation  error, 
or  roundoff  error  depending  on  the  scheme  used  to  obtain  the  fixed  binary 
number. 

For  two ' s complement  representation , the  error  Et  for  truncation  was  found  to  be 
0 2 Ex  > -2"^.  Rounding  off  a two's  complement  number  will  produce  an  error 
given  by  the  following  expression: 

-(l/2)(2-l>)  < Er  i (1/2)  (2-b) 

The  mean  and  variance  of  the  quantization  noise  are  shown  to  be 

■ 0 

Oe*  - A*/12  - 2“2b/i2 
for  rounding  and 
m^  - -2“b/2 

Og*  - 2-2b/i2 
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for  two's  complement  truncation. 

If  we  assume  that  the  signal  power  to  the  A/D  converter  Is  then  the 

signal  to  noise  ratio  produced  by  the  quantization  will  be  > 

(12/2-2b){j^i. 

Expressing  In  terms  of  db,  we  have  the  following: 

SNR  - 10  logjo  (Ox*/Oe*)  - 6b  + 10.8  + 10  logio  (0,^*) 

It  can  be  seen  that  the  signal  to  noise  ratio  Increases  6 db  for  each 
bit  added  to  the  register  length.  To  prevent  the  signal  from  exceeding  Che 
dynamic  range  of  the  A/D  converter,  Che  Incoming  signal  can  be  scaled.  Let 
A be  the  scaling  factor,  the  resulting  signal  to  noise  ratio  will  Chen  be: 

10  logio^^^2^y  6b  + 10.8  + 10  logio  (OjE*)  + 20  log^Q 

Since  A Is  less  Chan  1,  20  logj^Q  (A)  Is  negative  and  so  by  reducing  the 
Input  amplitude,  the  signal  to  noise  ratio  Is  reduced.  If  A Is  assumed  to  be 
(1/4)0,^,  the  signal  Co  noise  ratio  Is  (6b-1.24)db.  To  obtain  a signal  to  noise 
ratio  of  60  db,  10  bits  are  required. 

The  effect  of  quantization  for  Che  coefficients  Is  slightly  different 
from  that  of  the  A/D  conversion  process.  In  the  A/D  conversion  process,  noise 
Is  Introduced  Into  the  signal.  In  the  quantization  of  the  filter  coefficients, 
the  effect  Is  a change  from  Che  Ideal  frequency  response. 

A. 2 FIR  Filter  Noise 

The  filtering  process  of  a signal  Is  given  by  the  following  equation: 

N-1 

y(n)  ■ Z h(k)  X (n-k) 
k>0 

which  Is  a convolution  sum  process.  The  direct  realization  of  such  a system 
Is  a direct  realization  of  the  convolution  sum  process.  The  flow  graph  of  the 
realization  Is  shown  In  figure  A3. 

Since  the  multiplication  of  two  N bit  numbers  produces  a number  with  2N 
bits,  the  product  may  be  truncated  to  reduce  to  N bits.  This  corresponds  to 
the  addition  of  noise  to  Che  filter.  Assuming  that  each  product  Is  rounded 
off  before  summation,  the  output  noise  Is 

N-1 

f(n)  - Z eit(n) 
k-0 
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and  Che  variance  of  the  output  noise  Is 


Of* 


M 


2-2b 

12 


where  M ■ number  of  coefficients 

Notice  that  the  noise  Is  cumulative  and  depends  on  the  number  of  coefficients 

“!?  u*"*?  ^ computed.  Notice  also  that  If  all  2N  bits 
are  kept  for  addition  with  the  accumulator  output  truncated  to  N bit  the 

Che  A/D  conversion  process.  ^ 
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APPENDIX  B 


FIR  FILTER  DESIGN 


There  ere  many  ways  Chat  finite  lapulse  response  filters  can  be  designed. 
One  of  the  design  techniques  is  by  the  Fourier  Series  method  In  which  the 
Fourier  Series  representation  of  the  frequency  response  of  the  filter  Is  com- 
puted and  used  as  coefficients  for  the  filter.  Since  the  Fourier  series  have 
Infinite  representations,  these  coefficients  have  to  be  truncated.  However, 
direct  truncation  leads  to  the  Gibbs  phenomenon  which  manifests  Itself  as  a 
fixed  percentage  overshoot  and  ripple  before  and  after  an  approximated  dis- 
continuity In  the  frequency  response.  The  size  of  these  overshoots  does  not 
reduce  by  Increasing  the  Impulse  response  duration.  Instead,  the  overshoot 
la  confirmed  to  a smaller  and  smaller  frequency  range  as  N Is  Increased.  To 
solve  this  problem,  various  functions,  known  as  windows  are  used  to  smooth 
the  filter  response.  The  various  windows  are  shown  In  figure  Bl. 

While  the  window  design  technique  Is  good.  It  has  Its  flaws.  One  problem 
Is  that  one  needs  a closed  form  expression  for  the  Fourier  Series  coefficient 
which  may  be  complicated  to  obtain.  A more  Important  problem  Is  that  the  band 
edges  cannot  be  specified  exactly  since  the  window  "smears"  the  discontinuity 
In  frequency. 

A second  method  for  designing  FIR  filters  Is  known  as  the  frequency 
sampling  method.  In  this  method,  the  frequency  response  of  a filter  Is 
specified  and  samples  are  taken  at  equal  frequency  spaclngs.  The  Impulse 
response  Is  obtained  by  using  the  interpolation  equation.  In  the  frequency 
sampling  method.  In  order  to  obtain  the  best  result  It  is  Important  that  the 
value  for  the  transition  region  be  an  optimum  value.  To  do  the  optimization 
problem,  a series  of  equations  have  to  be  written  and  solved  mathematically 
for  the  solution.  This  Involves  linear  prograsmlng  and  a computer  Is  needed 
as  mil  as  a simplex  algorithm  for  minimization.  The  disadvantage  to  this 
design  technique  Is  that  the  selection  of  N as  well  as  the  selection  of 
transition  samples  are  based  on  experience  and  requires  Iterations  of  the 
computer  program  to  find  the  optimal  transition  samples.  Notice  also  that 
the  rlpplea  cannot  be  specified. 

A third,  and  recommended  technique  for  designing  digital  FIR  filters  is 
by  means  of  the  Remez  Exchange  Algorithm.  The  filters  designed  by  using  this 
technique  have  proved  to  be  optimal  filters.  Optimal  In  the  sense  that  the 
peak  approximation  error  over  the  entire  Interval  of  approximation  Is  mini- 
mized. This  program  was  run  on  the  NADC  CDC  6600  computer  and  used  for  the 
design  of  the  various  filters  and  differentiator  In  the  algorithms.  The 
program  Is  listed  In  the  following  pages. 
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=«OGi»»V  SlE^EHC  CDC  bbDO  FT^  -/3.0-P3SC  CPT  = 1 79/03/29.  13 .2  8.6b . 

SPCGRAN  RcbEXC  (INPUT, CUTPUTt 
C 

c 

c 

5 C FBCCRAP  FOR  The  OESIGN  OF  LINEAR  PHASE  FINITE  IHPULSE 

C RESPONSE  (FI»»  filter  USING  THE  RtFEZ  EXCHANGE  AlCORITHH 
C three  TYF£S  of  FILTERS  ARE  INCLUOED— EANCPASS  FILTERS 
C DIFFERENTIATORS,  ANO  HILOERT  TRANSFORH  FILTERS. 

C 

10  C THE  INPUT  DATA  CONSISTS  OF  6 CAROS. 

C 

C CARO  1--FILTER  LENGTH,  TYPE  OF  FILTER  (1-HULTIFLE 
C PASSBAHO/STOPBANO,  2-OIFFERENTIATOR.  3>HILBERT  TRANSFORM 

C FILTER!,  NUH9ER  OF  BANCS,  ANC  GRID 
15  C OENSITY.  (FORMAT  - bill 

C 

C CARO  2 — EANOECGES,  LONER  ANO  UPPER  EDGES  FOR  EACH  BAND 
C NITH  A MAXIMUM  OF  « BANOS.  (FORMAT  - IBFS.bl 

C 

20  C CARO  3— CESIREO  FUNCTICN  (OR  DESIREC  SLOPE  IF  A 

C CIFFERENTIATORI  FOR  EACH  BANO.  (FORMAT  - 8F5.2) 

C 

C CARO  b — HEIGHT  FUNCTION  IN  EACH  BAND.  FOR  A 

C DIFFERENTIATOR,  THE  HEIGHT  FUNCTION  IS  INVERSELY 

25  C PRCPORTlCNAt  TO  F.  (FORMAT  - 8F5.2) 

C 

C CARO  5 — RIPPLE  ANO  ATTENUATION  IN  FASSBANC  AND  STOPBAND. 

C THIS  CARO  IS  OPTIONAL  ANO  CAN  BE  USED  TO  SPECIFY  LOHPASS 

' C filters  DIRECTLY  IN  TERMS  OF  paSSBANO  RIPPLE  ANO  ST0P6ANC 

so;  C ATTENUATION  IN  OB.  THE  FILTER  LENGTH  IS  OcTERNINEO  FROM 

C THE  APPROXIMATION  RELATIONSHIPS  GIVEN  IN  t 

C L.  R.  RAEINER,  approximate  OESIGN  RELATIONSHIP  FOR 

C LOHPASS  CIGITAL  FILTERS,  IEEE  TRANS.  ON  AUDIO  AND 
C ELECTROACUSTICS.  VOL.  AU>21,  NO.  5,  OCTOBER  73. 

35  C ••••HHEN  THIS  OPTION  IS  USED  THE  FILTER  LENGTH  ON 

C CARO  1 SHOULD  BE  SET  TO  O.****  (FORMAT  - 2015.81 
C 

c sample  IKPUT  OATA  SETUO 

C 

bO  C the  FOLLOHING  INPUT  OATA  SPECIFIES  A LENGTH  32  8ANCPASS 

C FILTER  HITH  STOPBANDS  9 TO  9.1  ANO  0.b25  TO  0.5,  ANO 

C PASSBANO  from  9.2  TO  0.35  HITH  HEIGHTING  OF  10  IN  THE 

C STOPBANDS  ANO  1 IN  THE  PASSBANO. 

C THE  GRID  CENSITY  IS  32. 

b5  C 932001003032 

C .0000. 1000. 2000. 35CO.b250. 5000 

C .00001.000.9009 

C 10.001.00010.00 

C 

50  C The  FOLLOHING  INPUT  OA^A  SPECIFIES  4 LENGTH  32  HlDcBANO 

C DIFFERENTIATOR  HITH  SLOPE  1 AND  HEIGHTING  OF  1/F.  THE  GRID 

C OENSITY  IS  ASSUMED  TO  BE  16. 

C 032002001000 

C .0000.5009 

55  C l.COO 
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C 1.30C 

C 

C TMc  FOULOWING  input  CA^a  specifies  a LOWPASS  filter 
C at  GIVING  CIRECTLT  TmE  RIFFLE  AND  ATTENUATION. 

C 901001002032 

C .0000. 20  00. 3009. SOCO 

C l.OOO.COOO 

C l.COOl.OCO 

C C. 200000000000940. COOOCOOCOuOG 

C 

CCNNCN  OES.MT, ALPHA,IFyT,NFCNS.NGi;iC«PI2.AO.CEVtX,r  .GRID 
CIH£NSI0N  IEXTIEG) ,A0(e6itALFHA(66>tX(66ltY(66> 

OIHENSION  HIM  i 

ai FENS ION  OESI 1045), GRID (1045) , NT (1045) 

OIHENSION  EDGE (20) . fx ( 10 ) , MTX ( 10 ) .OEVIATIIO ) 

CIHENSION  OHEGA(50) fRESPAISO  ) 

OOUSLE  PRECISION  ONEG A , SUN AR , SUHAC.RESP A , ATTN , R IPPlE ,SVRIP, CELT AF 
DOUBLE  PRECISION  PI2.PI 
COUBLE  PRECISION  AO.OEV.X.Y 
FI2=6. 28318530 7179586 
FI=3. 1415926535899793 
C 

C THE  FROGRAF  IS  SET  UP  FOR  A HAXIHUH  LENGTH  OF  128,  BUT 

C THIS  UPPER  LIMIT  CAN  BE  CHANGED  3Y  REOIHENSIONING  THE 

C ARRAYS  lEXT,  AO,  ALPHA,  V,  Y,  H TO  BE  NFMAX/2  * 2. 

C the  arrays  OES,  GRIC,  and  MT  must  6£  DIHENSIONEO 
C 16(NFHAX/2  42). 

C 

NFFAXsl2e 
190  CONTINUE 
JTYPE=0 
C 

C PRCGRAP  INPUT  SECTION 

C 

REAOl.NFILT, JTYPE.NBANCS.LGRIO 

1 FORMAT  (4131 

lF(NFlLT.GT.NFHAX.OR.NFILT.Ea.3)  CALL  ERRCR 
IFINEANOE.LE.O)  NPAN0S=1 

c 

C GRID  CENSITY  IS  ASSUMEC  TO  BE  16  UNLESS  SPECIFIEO 

C OThERNISE 

C 

IFILGPIO.LE.OI  LCRI0=1E 
je»2*NBANCS 

REA02,  (EOGE(J) ,J=1,16I 

2 FORMAT  (16(F5.4)) 

READS,  IFX(Jt , J3i, ai 

3 FORMAT  (e(F5.2)) 

REA04,  (HTX(J) , J«l, 8) 

4 FORMAT  1 8 (F5.2)  ) 

IFIJTYPE.EO.O)  CALL  ERROR 
IFINFILT .EQ,Q)  GO  TO  19 
NROX>C 

GO  TC  11 

10  REA05, RIPPLE. ATTN 
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PROGRAM 

115 

120 

125 

130 

135 

lUO 

145^ 

150 

155 

160 

165 


"Et'EXC  COC  6600  FTN  V3.0-P380  CPT  = l 79/03/29.  13.28.54.  FAGc 

5 FORMAT  (2015.81 

riffle*! '’••IRIFOLE/20) 

ATTN»1C**(-ATTN/20) 
aiFPLE*«RIPPLE-l)/ (RIPFLt+l> 

RIFPLE*OLCG10 (RIPPLE) 

ATTN*CL061C(ATTN) 

SVRIPsRIFFcE 

R IPPL£  = ( 0. C053C 9»(RIPPLE**2)»C. 0711 4*RIPPLE-C  .4761) •ATTN 
l*(-0.00266*(aiPPLE**2) -0 .5941*RIPPLe-0 . 4278) 

ATTNsll. 01217*C.51244*SVRIP-0.51244*ATTN 
CELTAF*ECGE(3I-EDGE(2) 

NFILT=fiIFFLE/DELTAF-ATTM*(0eLTAF**2)*l 
NROXsl 
11  MEG=1 

IF(JTYPE.EQ.l)  NEG=0 
N000*NFlLT/2 
NOCO=NFILT-2*NOOO 
NFCNSsNFILT/2 

IF(NOOC. EO.I.ANO.NcG.EQ.O)  NFCNSsNFCNS^l 
C 

C SET  UP  THE  GRIO  DENSITY.  THE  NUMBER  OF  POINTS  IN  THE  GRID 
C IS  (FILTER  length  ♦ D’GRIO  OfNSITY/2 
C 

GRIOdlsEOGEd  > 

OELF=LGRID*NFCNS 
CELF=C.5/CELF 
IFtNEG.EQ.u)  GO  TO  135 
IF(EDG£(1).LT.0ELF)  6RI0(1)*0£LF 
135  CONTINUE 
J*1 
L=1 

LBAND*! 

140  FUFsECGE (L^l) 

145  TEMPsGRICIJ) 

C 

C CALCULATE  THE  DESIRED  MAGNITUDE  RESPONSE  AND  THE  HEIGHT 
C FUNCTION  ON  THE  GRIO 
C 

CES(J)=£FFCTEMF,FX,NTX»LEAN0. JTYO£> 

HT(JI*HAT£(TEHP,FX tNTX.LEANOtJTYPE) 

J=J*1 

GRI0(JI*TEHP«0ELF 
IF(CRI0(4) .GT.FUP)  GO  TO  150 
GO  TC  145 
150  GRIO)J-l)*FUF 

0£S(J-l)*EFF(FUP.FX,HTX.L8AN0t JTYPE) 
NT(J-l)>HATE(FUPtFXtHrXt LBANOt JTYPE) 

LBAN0*LBANC^1 

L»L^2 

IF(LBANO.GT.NBANOS)  GO  TO  160 
CR10(J)*£CG£(L) 

GO  TC  140 
165  NGPIC=J-1 

IF(NE6.N£.N000)  GO  TO  165 

IF(GRIO(NGRID) .GT . (0 .5 >0£LF) ) NCR 10*NGR lO-l 
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170 


175 


100 


185 


19C 


195 


200 


205 


210 


215 


220 


16*  CONTINUE 

C 

C SET  UP  A NtH  APPROXIMATION  PROOLEM  WHICH  IS  RQUIWAlENT 

C TO  TH&  OftlGlNAk  PROBLEH*  . 

C 

IFINEC)  170,170,180 
170  IFtNCOO. EQ.1>  GO  TO  200 
00  175  J^i.NGPIO 
CHANGE =0 COS (PI *GRI0(J> i 
. OES(JlsOeS(J>/CHANGE 
175  WTIJIawr (JI*CHANGE 
GO  TO  200 

len  IFINCCC.EO.il  GO  TO  190 
00  185  J=1,NGRIC 
CHANGE<OSlN(Pt *0910(31 1 
0ES(3)*0:S(J>/CHANCE 
10*  HT(J|3HT (J)*ChANGE 
GOTO  200 

190  CO  195  J*1,NGR1C 

CHANGEp0S1N(PI2*GPI0( Jl I 
0£S(J)pOES(J)/CHANG£ 

195  WT(JtsHT(J)*CHANGE 
C 

C INITIAL  GUESS  FOR  THE  EXTREMAL  FREGUENCIES— EQUALLY 
C SPACED  ALCNG  THE  GRIO 
C 

200  TEHP*FtOAT|NGRIO-ll/FLCAT(NFCNS) 

CO  210  Jsl.NFCNS 
210  lEXT(J)* (J-1I*TEHP*1 
IEXT(NFCNS*l|sNGRIO 
NMlaNFCNS-l 
NZsNFCNS*! 

C 

C CALL  THE  PEME2  EXCHANGE  ALGORITHM  TC  00  THE  APPROXIMATION 
C PROBLEM 

C 

CALL  RENEZ(EOC£,NBANOS) 

C 

C CALCULATE  THE  IMPULSE  pESPONSE. 

C 

IF(NEG)  300,300,320 
300  IFINOCO.EO.OI  GO  TO  310 
00  305  J*1,NH1 
305  H(J|s0.5*AlPHA(NZ-JI 
• H(NFCNS>*ALPHA(1) 

GO  TO  350 

310  h(l)s0.2**ALPHA(NFCNSI 
00  315  J32,nhi 

315  HI3I>0.25*(ALPHA(NZ-JI  *ALPHA(NFCNS*2'>JI  I 

HINFCNSI =0.5*ALFHA (1)*C. 25*ALPHA(2» 

GO  TC  351’ 

320  IFINODO.EC.O)  GO  TO  330 
HI1)«0.25«ALPHA(NFCNS» 

H(2I«0.25»AlPHA  (NMD 
00  325  Js3,NMl 
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325  M« j»a0 .25*(ALPhAJNZ-J» -ALPHA (NFCNS^3-j» ) 

HINFCNSI =0.5*ALPHA (1) -r,25*ALPMA( 3) 

HINZt«C.O 
GO  TO  350 

225  330  H<1)30.25«ALPHA(NFCNS) 

00  335  J:2.NN1 

335  HCJ|cg.25*IALPHAINZ-J)-ALPHA(NFCNS«2-JI I 
H(HFCNSI 3C.5*ALPHA<1I-0.25«ALPHA(2) 

C 

230  C F9CGAAP  OUTPUT  SECTION 

C 

PRINT  3V0 
3b0  FORNAT(lHl) 

350  PRINT  36" 

235  360  FOPHAT(25X, •FINITE  INFULSt  RESPONSE  (FIRI*/ 

125X, •LINEAR  PHASE  DIGITAL  FILTER  DESIGN^/ 

225X.^RENEZ  EXCHANGE  ALGORITHM^/) 

IF(JTYPE.EQ.1»PRINT  365 
365  FORHAT(25X, •BANDPASS  FILTER^/) 

240  IF  (JTVPE.EQ.2)  PRINT  T’O 

37P  F0PHATI25X,^0IFFER£ntIAT0R*/> 

IF(JTYPE.£Q.3)  PRINT  375 
375  FORMATC25X, •HILBERT  TR ANSFORNER^/ I 
PRINT  37«»NFILT 

245  378  FORHAT<25X, •filter  LENGTH  a •,I3/» 

IF(NROX.EO.l)  PRINT  3'»o 

375  FOfiHAT«25X,^FILTEP  LENGTH  OETERHINEO  BY  APPROXIHATION^/) 
PRINT  383 

380  FOPHAT(23X,^ INPULSE  RESPONSE  •» 

250  00  381  Jal.NFCNS 

K=NFILT*1-J 

IF(NEG.EQ.CI  PRINT  3S2,J,H|J),K 
IF(NEG.EC.l)  PRINT  303«J,HIJ),K 

381  CONTINUE 

255  382  FORMAT(20X,^HI^,I3,^I  * •,615. 8. • * H(^,I4,^)^» 

383  F0RHAT«2CX,^HI^,I3,^I  » •,E15.8,^  * -H<^,I4,^I^» 

IFINEG. EC. 1. AND. NOOO.EO.il  PRINT  384, NZ 

384  FORMATJ2CX,^HI^,I3,^>  * O.O^I 
00  450  Kal,NaAN0S,4 

260  KUPsK«3 

IFtKUP.GT.NBANCSI  KUPaNBANOS 
PRINT  385,(J,JaK,  KUPI 

385  F0RNAT(/24X,4l^aAN0^,I3,8XI) 

PRINT  390, (EaG£(2^J-l) ,JaK, KUPI 

265  390  FORHATI2X, •LONER  9AN0  iOCE^, 5F19. 91 

PRINT  395, IE0CE(2^Ji,JaK, KUPI 
395  FORMAT (2X,^UPPER  BAND  EOGE^, 5F19. 91 

IF(JTYPE.NE.2)  PRINT  400 , (FX  IJi . J«K,KUP) 

400  FCRHAT(2X,^0ESIRE0  VALUED, 2X,5F15, 9) 

270  IF<JTYP£.E0.2)  PRINT  4"5 , < FX ( Jl  , J»K,KUF I 

405  FORHAT«2X,*OESIREO  SLOFE^, 2X , 5F15 . 5) 

PRINT  413,(HTX(JI,J«K,KUP) 

410  F0RHAT(2X,^NEIGHTING*,6X,5F15.9I 
CO  420  JsK,KUP 

275  420  OEVIATUIaOEV/HTXIJI 
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260 


26! 


290 


29S 


300 


305 


310 


PRINT  <i25,«05VIAT(JI,J  = K.KUP» 

425  FORHAT »2X, •DEVIATION*, 6X,5F15.9» 

IP<JTYPE  .N€.U  GO  TO  450 
00  430  J:K,KUP 

IF  (FX IJ)  .EC.l.Ot  OEVIAT(JI«(1.0*OEVIAT(Jli/(t.Q>OEVIAT(JII 
430  OcVIATiJI s20.0*AlOG19(OevlAT(J» ) 

PRINT  435,IDEV1ATU),JsK,KUP) 

415  FOfiNAT(2X, •DEVIATION  IN  08*,5F15.9» 

450  CONTINUE 

PRINT  455, (GRICtlEXTC Jl) .Jsl.NZ) 

455  FOPHAT</2X,*£XTRENAC  FREOUENCI £$•/ e2X , 5 F12. 7 1 ) 

PRINT  460 
460  FORNATIlNll 

C CALCULATE  FRcQUENC/  RESPONSE 

PRINT  710 

710  FOPHATUax* FREOUENCT  RESPONSE  •> 

CO  610  IKA  s 1,50 

CNEGA(IKA)  s 0 .0100000GO*I1KA  - II 

SL^AR  « 1,0 

SU^iAC  X O.C 

CO  620  NIK  S 1,  NFCNS 

SUNARxSUHARaHINIKI •C0CCS(PI2*0HECAIIXAI  •<NIK-1) I * 
10COS(PI2«ONEGA IIKA) •INFIlT-NIKI  1 1 
SUMAC  s SUMAC  • H(NIK)*(0SIN(PI2^0NEGA(IKA)*(NIK-1)I  t 
10SINIPI2*ONEGA(IKA|*(NPlLT-NIKni 
620  CONTINUE 

IF  INCOO.EO.ll  SUNAPaSUNAR-HINFCNSI'OCOSIPIExCNEGAlIKAI^INFCNS 
1«1>) 

IF  (NOOO.EO.l)  SUHACxSUNAC>HINFCNS>*OSlNIPI2^CNEGA(IKAI*fNFCNS 
I'll  I 

RESPA (IKA)  « 0SaRT(SUNAR**2  • SUNAC**2) 

RESPAdKAI  s 20.0^0LOG10(RESPA(IKAI) 

PRINT  630,  ONEGAdKAI  , RESPAdKAI 
630  FORMAT (15X, 012. 4, 1SX,C12.4} 

610  CONTINUE 

IFINFILT.NE.O)  GO  TO  130 

CONTINUE 

ENC 
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TION  £®F  COC  66uC  FTN  V3.D-P3aC  CPT*l  79/>)3/23.  13 .2  d.  64, , F4CE 

FUNCTION  £FFTTEmb,FX,HTX,LBANO, JTVFEJ 

FUNCTION  TO  CALCULATE  THt  OESIREO  HACMTUOE  RtSFONSE 
AS  A FUNCTION  OF  FREQUiNCT 

OIHENSION  FX<SI,HTXI9I 
IF(JTVP€.E0.2|  CO  TO  1 
SFF»FX(LeANO) 

RETURX 

1 £FF«FX CLeANO)»T£HP 
RETURN 
ENO 
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H»TE  COC  6600  FTM  V3.0>PS6P  OPT»i  79/03/29.  13 .2 8 .5<. . p«CE 
FUNCTION  NATsJTENP.FX.NTX.LBANO.JT  fF£» 

FUNCTION  TO  CALCULATE  THE  NEIGMT  FUNCTION  AS  A FUNCTION 
OF  FREQUENCY 
OIHENSION  FX(9t,HTX(5) 

IF(JTVPE.EC.2)  GO  TO  I 
NATEshTX(LBANO> 

PETUPN 

1 if(fx«lbanoi.lt.o.oci)  go  to  2 

MATEaNTX (LdANO) /TEHF 
RETUNN 

2 MATE«NTX  (LBANOI 
RETURN 

ENO 
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REI>£Z  COC  66CC  FTK  yS.H-Pita  OPT*l 


79/C3/29.  13. 2S.3<..  PACE 


1 


5 


10 


15 


20 


25 


30 


35 


<tO 


A5 


50 


55 


subroutine  REHEZIEOCE.NBANOSI 

r 

c this  subroutine  INPLENENTS  The  RcHEZ  exchange  algorithm 

C ^OR  the  McIGHTEO  chebychev  approxihation  of  a continuous 

c function  with  a sun  of  cosines,  inputs  to  the  subroutine 
C are  a cense  grid  MHICH  replaces  the  FRcAUENCV  axis,  the 

C OESIREO  FUNCTION  ON  THIS  GRID.  THE  HEIGHT  FUNCTION  ON  THE 

C GRID.  THE  KUNBER  OF  COSINESt  ANC  AN  INITIAL  GUcSS  OF  THE 

C EXTREMAL  FREQUENCIES.  THE  RROURAN  NININIZES  THE  CHE3YCHEV 

C ERROR  BY  CETERHINING  THE  BEST  LOCATION  OF  THE  EXTREMAL 

C FREQUENCIES  (POINTS  OF  HAXINUN  ERROR*  ANC  THcN  CALCULATES 

C THE  COEFFICIENTS  OF  THE  BEST  APPROXIHATION. 

C 

COMMON  OEStHT.ALPHA.IEXT.NFCNS.NGRIO.PIZ.ACtOEV.X.Y.GRlO 
OIMENSIOh  EDGE (23* 

OIHEASION  lEXTIEbl « AOI 6GI • ALPHA (66 t «X (66*.Y(66> 

OIMENSION  OESdCNSI  tGRIO(lCNS)  .NT(1569l 
OIHENSION  A(66ltP(65ltQI65* 

OOUSLE  PRECISICN  PI2,0NUM« OOEN.OTENP* A, P, Q 
COUBlE  precision  AC.OEV.X, Y 
C 

C THE  PROGRAM  ALLOWS  A HAXiNUN  NUMBER  OF  ITERATIONS  OF  25 
C 

ITRNAX>25 
CEWL»-l.lJ 
NZsNFCNS^l 
NZZ«NFCNS«2 
NITER»C 
100  CONTINUE 

1EXT(NZZI>NGRIC«’1 

NITER«NITER»1 

IF(NITER.GT.ITRHAX)  GO  TO  AuO 
CO  lie  Jsl.NZ 
OTEMPsGRIOdEXTCJII 
OTEMF»OCOS(OTEMP*PI2> 

110  X(J»*OTEMP 

JETm(NFCNS>1)/15»1 
CO  120  J«ltNZ 
120  AC(JI>0( J.NZ.JcT) 

CNUNsO.O 

CC£N«0.0 

K«1 

00  130  Jsl.NZ 
L«IEXT(J) 

CTEMF>A0(JI*0ES(LI 

CNUn>CNUM«OTEMP 

CTEHP«K*AO(JI/MT(l) 

OC£N«COEN«OTEMP 
130  K«>K 

CEW«CNUM/CCEN 

NUsl 

IF(OEV.CT.i}.Ot  NUs>l 

CEO«-NU*C£Y 

K«NU 

CO  160  JalfNZ 


I 


B-12 


NADC-791S3-40 


SUBROUTINE  a£l>E2  COC  6E0C  FTN  V3.0*Pi80  0PT>1 


79/03/29.  13. 24. SR. 


60 


65 


70 


75 


SO 


05 


90 


95 


100 


1C5 


110 


. L»r£XTIJI  P4C£ 

CT£NF»IC*0EW/NT«L> 

/fJI ■0£S(Lt  »CT£nP 
IbO  K«*K 

IF(0£V.GE.3EVLt  GO  TO  ISS 
CALL  CUCH 
GO  TO  i.Q)' 

150  C£tfL«CEV 
JCHNGEsO 
KlsIEXTIl) 

KN2«IEXTINZI 

KLOMmO 

NUTs-NU 

J*1 

C 

C SEARCH  7CA  THE  EXTREMAL  FREQUcNClES  OF  THE  BEST 

C APPROXINATION 

C 

200  IFU.EQ.N2ZI  YN2«COMP 
IF(J.GE.h2Z>  GO  TO  300 
XUP«XEXT  U*l> 

LsIcXTtJI*! 

NUT--NUT 

IF(J.EQ.2I  YlsCQMP 
COMP>CEV 

IFIL.GE.KUP)  GO  TO  220 
ERR>GEE(L,NZt 
ERR>(ERR>OESILn*NT(Lt 
OT£NPsNUT«ERR>CONP 
IFlOTEHP.LE.O.Ot  GO  TO  220 
COHP«NUT*ERIt 
210  LsLRl 

IFlL.GE.KUPI  GO  TO  215 
£RR«GEEIL«N2) 

£RR«(ERR-0£S(LI IPHTIL) 

OTENP«NUT»£RR-COMP 
IFtOTENP.LE.O.OI  60  TO  215 
COMP* NUT ’ERR 
GO  TC  21P 
215  IEXTCUML‘1 
J«J*1 
KLCM*L-1 
JCHNGEsJCHNGlR 1 
GO  TC  200 
220  L*L-1 
225  LaL‘1 

IFIL.LE.RLONI  CO  TO  25'* 

ERR«GE£(L«N2t 
ERR>(ERR-OcSILi >«NT(l) 

CT£MP«NUT*£RR-CONP 
IF(0TEHP.CT.0.C»  GO  TO  230 
IFCJCHNGE.LE.Oi  GO  TO  225 
GO  TO  260 
230  COMpaNUT*£RR 
235  L»L-1 


2 
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SUBBOUTTNE  «£H£2  COC  66(IC  FTN  V3.0-P380  0PT«1 


79/03/25. 


13. 28. $4.. 


115 


120 


125 


130 


135 


-0^0 


1<.5 


ISO 


155 


160 


165 


IF(C.Lt.HLCWI  GO  TO  24!’ 

£RP>GE£(L,NZ) 

Ei;R«IERR«OESail«HT(ui 

QT£Mf«m»T*£RR-COHP 
irii)TENP.LE.O.«>  GO  TO  24C 
C0NP>NU14EPR  * 

GO  TC  235 
240  KLO»(»I£X  T«j| 

IExT(ji«l»i 

JS  J»1 

JCHRGE*3CHNG£« 1 
GO  TO  200 
250  L«l£ITIJI«l 

IF<jChH6i.G’.&>  CO  TO  215 
755  L»t*l 

IFIL.CE.KUPI  CO  70  260 
E«R«GEEUtN2) 
£«li>l£ltR>OcS(CII»MT(LI 
0 T s HP»  NUT*ERR-CONP 

If iotekp.le.o. ai  co  to  255 
COPP«»*UT»tRR 
GO  TO  210 
26C  KLCN«1EXTCJ> 

J«  J4l 
CO  TO  200 

300  IflJ.CT.42ZI  CO  TO  320 

IFtKl.GT  .lEXTI 11)  Kl«r7XT(l» 

ZFIKNZ.LT.ieXTlNZl I KNZ*tEXTINZ) 

NUTlsMUT 

NUTb'NU 

L*" 

KUP*X1 

CONP»TN2»«1.POCOU 

tUCK*! 

310  L*L»1 

IflL.CE.KUPI  GO  TO  315 
ERI(>GECU,NZ) 

:RR«<ERR>OES(CI I’RTTUI 
3TcHP*NU T»£RR-COHP 
IF  IOTENP.l.e.3.  Ci  CO  TO  310 
COPP«RUT»£PR 
J*NZZ 
CO  TO  210 
315  LUCK«6 

GO  TC  325 

320  IFILICK.GT.9I  CO  TO  357 
IFICCPP.CT.Yl)  Tl«COHP 
KlsIEXTIRZZI 
325  L>NGRIO«l 
KLOMzKNZ 
NUT*-mjTl 
COPPaYl* (5.030011 
330  L«l-1 

IF(L.L£.KL0M»  CO  TO  345 
£RP<GEE(L,NZ) 
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o*H€Z  COC  abfci'  FTh  W3.0-P38P  0PT»1  79/0  3/29.  1 3.2S . S<>  .P«G£ 


1T« 


175 


180 


185 


190 


199 


200 


205 


210 


215 


220 


'KS»(EBR-DES(L»I»MT(L» 

CTEHF»NUT»i«9-C0HP  ’ 

IFIOTEMT  .Lc.O.OI  GO  TO  330 
jxK22 

COPP»*tOT»ERI» 

LUCKS*.  UCK»10 
CO  TO  235 

’SfcO  IFILlCK.EO.ei  GO  ’0  T’O 
* CO  3«.5  Jsl.NFCNS 

31.5  IEXTINZ2>J|sI£X7(NZ-JI  . 
lEXTIltsKl 
GO  TO  10" 

S«T  KNsieXTINZZt 

CO  360  Jsl.NFCNS 
360  I-XTtJI»r*XTU*ll 
lEXTIKZ) iKN 
GO  TO  10" 

370  IFUCHNCE.CT.ai  GO  TO  ItO 

C CAICULATICN  OF  TMR  COEFFICIENTS  OF  TNE  BEST  APPROXIMATION 
C USING  THE  inverse  DISCRETE  FCURIER  TRANSFORM 

1.00  continue 

NH1«NFCNS«1 

FSM-l.OE-06 

GTFMPsGRICTlI 

X|NZZ»»-2.? 

CN»2»NFCNS-1 

CELF«1.0/CN 

L»1 

KKK»e 

IFtEOGEIl).EQ.0.0.ANa.ECG4l2«NaANOSi.EO.O.5l  KKKsi 
IFINFCNS.I.E.31  KKKsl 
IFIKKK.EC.1)  GO  TO  609 
CTENP«0CCS(PI2*GRICI1»» 

ONUMsOCOS (P12»ORIO INGRID) ) 

AA>2.0/(ST£HP-CNUM| 

BB«-IOTEMP*ONUM>/(OT-HP-ONUN) 

NCR  CONTINUE 

CO  630  JslfNFCNS 

FT*(j-n*C£LF 

XT*0CCSCP12*FT) 

IFtKKK.EO.l)  CO  TO  610 
XTs|XT*BBI/AA 
FT*ACCSIXTI/PI2 
61C  XEsXIL) 

IFIXT.GT.XEI  GO  TO  620 
IF((Xi-XT).t.T.FSH)  CO  TO  615 
L»L»1 
GO  TC  619 
*15  AUI*TIL> 

CO  TC  62E 

*29  IFnXT-Xi).LT.FSH)  GO  "0  615 
GRlOllIsFT 
A<J) «C£E(ltNZ) 

625  CONTINUE 

IFCC.GT.l)  L*C-1 
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subroutine  REf'EZ  COC  6fcOC  F TN  W3.J>P3aO  OPTal  7^/i3/29.  13. 24. So. 

<»30  CONTINU£  5 

G«I0<1>»GTS’HP 
CCiNxPI2/CN 
CO  Sic  Jxl.NFCNS 

225  0TEMP*C.3 

CNUM»IJ-H*03EK 
IFINHl.LT.il  CC  TO  EC5 
CO  5CC  K*1.NN1 

Sir  CTeNFsCTeHF*ACK*ll ‘CCOS (ONUH*Kl 

230  305  0T£HP«2.5*0T£MF»*«1I 

510  ALPHA! J> <CTENP 
OC  550  JX2.NFCNS 
550  ALFHA(J132*ALPHAIJ)/CN 
ALPHA! II : ALPHA! It /CN 

235  IF!KAK.EC.1)  go  TO  51.5 

P ! 11 *2 .0 • ALPHA !NFCNSI * cB  xALPHA  !NH1 1 
P!2l>2.a« AA*ALFHA!NFCN5) 
a!l)«ALPHA!NFCNS'2»>ALPHAINFCNSI 
CO  5N0  J«2.NH1 

2A0  IFIj.LT.NPil  CC  TO  515 

AA*0.5*AA 
3BsC.S*Be 
515  CONTINUE 
F!J'»1IeO.O 

2N9  CO  520  Kxl.J 

A!K|sFIKt 

520  P!K| x2.0*Ee*A!KI 

P!2)*F!2l*A!ll*2.3*AA 

JHlsJ-1 

250  00  525  Kxl.JHl 

525  P!K|aP!Ki«aiKI «AA*A!K*2I 
jPlaJ*! 

CO  5!C  KX3.JP1 
53:  PIKI«F!K)*AA*AIK-1I 

255  IF!J.EC.NH1I  GO  TO  5<.5 

00  5?5  K<1.J 
535  C!KI«-AIK) 

C!1|zQ!i I«ALPHA<HFCNS*1>JI 
5i>0  CONTINUE 

260  CO  5fc3  Jal.NFCNS 

5N3  ALPHA!JI«P!J1 

51.5  CONTINUE 

IFiNFCNS.GT.3l  RETURN 
ALPHAINFCNSTII  «0.0 

265  ALPHA!NFCNS*2I >0.0 

return 

ENC 
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FUNCTION  C COC  660C  FTN  V3.0-P380  CPT=1  7S/ta3/29.  13.2S.5n.  paG£ 


COUBLE  PRECISION  FUNCTION  OIK.N.NI 
C FUNCTION  TO  CALCULATE  THE  LAGRANGE  INTERPCLATICN 
C COEFFICIENT  FOR  USE  IN  THE  FUNCTION  GEE 

CONHCN  OES.HTt ALPHA,IEXTfNFCNStNGR10*PI2 tAOtOFV«X»V,CRIO 
S OlHENSION  IEXT(66>.A0(66l«ALPHA<66)>X<66)tYC66) 

CIHENSION  OES(1046),GRIO(10N5>  <HT(1045) 

COUBLE  PRECISION  AG.OEVtX.V 
CCUBLE  PRECISION  Q 
COUBLE  PRECISION  PIE 

10  031.0 

QsXCK) 

00  3 L31.M 
CO  ? JsLtN.H 
IF«J-KJ1,2,1 

15  1 Os2.0*D*  (u-xun 

2 CONTINUE 

3 CONTINUE 
0»l.C/0 
RETURN 

20  END 
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PUNCTICN  G??  COC  &6C5  FTH  V3.0-P580  OfT*l  79/<)3/29.  13. 26. PAGE 

COUBLE  PRECISION  FUNCTION  CE£(<.NI 
C FUNCTION  ’O  EVACUATE  TMt  FREQUENCY  RESPONSE  USING  THE 

C LAGRANGE  INTERPOLATION  FQRNULA  IN  THE  BARYCENTfilC  FORM 

CONHCN  OES.NT, ALPMA,IEXT,NFCNS.N6RlCtFI2.A0,0EW,X,Y,GRID 
5 OIPENSION  IEXT(66I «A0(66l«ALPHA(66l«X(66i«V(66l 

OINENSION  OESI10N9t .GRID (1 0451 .H^CIONEI 
GOueCE  PRECISION  P.CtO.XF 
COUBLt  PRECISION  PIZ 
COUBCE  PRECISION  AO.OEVtX.V 
10  PsO.O 

XFzGRIOf K) 

XF=0CCSJPI2*XF) 

0*0.0 

CO  1 J*1,N 

15  C=XF-X«JI 

C*AO(J)/C 
C*0«’C 

1 P=F*C*Y(J) 

GEE=F/D 

2(?  RETURN 

END 
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SUePOUTlNE  OUCH  COC  6600  FTK  W3.0-P380  OPT  = l 79/tt3/29.  13.28.5**.  f*CE 

subroutine  ouch 

fRlNT  1 

1 FOPHATI* failure  to  CONVERGE  

I’OFROBABLE  CAUSE  IS  MACHINE  R0UN0IK6  ERROR*/ 

2*0THE  IMPULSE  RESPONSE  HAV  BE  CORRECT*/ 

3*0CH£CK  NITH  A FREQUENCY  RESFONSE*) 

RETURN 

ENC 
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REPORT  NO 


NADC-79153-A0 


AIRTASK  NO.  WF21-200-000 
Work  Unit  No.  ZC902 


No.  of  Copies 

NAVAIRSYSCOM,  AIR-950D  5 

(2  for  retention) 

(1  for  AIR-370) 

(1  for  AIR-5330) 

(1  for  AIR-53301) 

CNO 2 

(1  for  OP-951) 

(1  for  OP-952) 

COMOPTEVFOR 1 

DDC  


